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ABSTRACT 

Motivated by recent spectroscopic evidence for carbon-rich atmospheres on some transiting exo- 
plancts, we investigate the influence of the C/O ratio on the chemistry, composition, and spectra of 
extrasolar giant planets both from a thermochemical-equilibrium perspective and from consideration 
of disequilibrium processes like photochemistry and transport-induced quenching. We find that al- 
O [ though CO is predicted to be a major atmospheric constituent on hot Jupiters for all C/O ratios, 

other oxygen-bearing molecules like H2O and CO2 are much more abundant when C/O < 1, whereas 
CH4, HCN, and C2H2 gain significantly in abundance when C/O > 1. Other notable species like N2 
and NH3 that do not contain carbon or oxygen are relatively unaffected by the C/O ratio. Disequi- 
librium processes tend to enhance the abundance of CH4, NH3, HCN, and C2H2 over a wide range 
of C/O ratios. We compare the results of our models with secondary-eclipse photometric data from 
the Spitzer Space Telescope and conclude that (1) disequilibrium models with C/O ~ 1 are consistent 
with spectra of WASP-12b, XO-lb, and CoRoT-2b, confirming the possible carbon-rich nature of 
these planets, (2) spectra from HD 189733b are consistent with C/O < 1, but as the assumed metal- 
licity is increased above solar, the required C/O ratio must increase toward 1 to prevent too much 
H2O absorption, (3) species like HCN can have a significant influence on spectral behavior in the 3.6 
and 8.0 /im Spitzer channels, potentially providing even more opacity than CH4 when C/O > 1, and 
(4) the very high CO2 abundance inferred for HD 189733b from near-infrared observations cannot 
be explained through equilibrium or disequilibrium chemistry in a hydrogen-dominated atmosphere. 
We discuss possible formation mechanisms for carbon-rich hot Jupiters, including scenarios in which 
the accretion of CO-rich, H20-poor gas dominates the atmospheric envelope, and scenarios in which 
the planets accrete carbon-rich solids while migrating through disk regions inward of the snow line. 
The C/O ratio and bulk atmospheric metallicity provide important clues regarding the formation and 
evolution of the giant planets. 

Subject headings: planetary systems — planets and satellites: atmospheres — planets and satellites: 
composition — planets and satellites: individual (HD 189733b, WASP-12b, XO-lb, 
CoRoT-2b) — stars: individual (HD 189733, WASP-12, XO-1, CoRoT-2) 



1. INTRODUCTION 

Transit and eclipse observations provide a unique 
means with which to characterize the physical and 
chemical prop erties of "hot Jupiter" e xoplanet atmo- 
spheres (e.g., iSeager and Demiirgl I2010I ). Of particu- 
lar importance in terms of furthering our understand- 
ing of planetary formation and evolution are the con- 
straints on atmospheric composition provided by multi- 
wavelength photometric and spectral observations of 
the star-planet system as the transiting planet passes 
in front of and behind its host star as seen from an 
observer's perspective. Transit and eclipse observa- 
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tions to date have been used to infer the presence of 
the molecules H2O, CO, CO2, and CH4 in the tropo- 
spheres and/or stratospheres of extrasolar giant plan- 
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Despite this remarkable feat, several factors can com- 
plicate the compositional analyses from transit and 
eclipse observations, and the results regarding molec- 



2 



Moses et al. 



ular abundances are typically model dependent. For 
example, in addition to being sensitive to the atmo- 
spheric composition, the transit and eclipse signatures 
are strongly affected by the atmospheric thermal struc- 
ture, the presence of clouds and hazes, the efficiency 
of atmospheric dynamics, and the atmospheric mean 
mole cular mass, bulk metallic ity, and elemental ratios 
(e.g., iSeager and Pemine] 120101 ) . The complexity of the 
problem is best attacked through a combination of so- 
phisticated t heoretical modeling and careful data analy- 
sis (e.g., seelMarlev et al.ll2007t iBurrows fc Ortonll2010l: 
Showman et all 120101 : ISeager and DemingJ l2010t iWinn I 



20101 and references therein), but degeneracies between 
solutions, and resulting uncertainties in composition, 
will continue to plague the determination of constituent 
abundances until higher-spectral-resolution observations 
become available to more unambiguously identify atmo- 
spheric species. 

The largest source of degeneracy arises from uncertain- 
ties in the atmospheric thermal structure and compo- 
sition. To get around this problem, one can attempt 
to derive simultaneous constraints on both the compo- 
sition and thermal structure through statistical treat- 
ments that systematically compare millions of forward 
models with the observed phot ometric and spectral data 
(|Madhusudhan fe Seagerll2009ft or through parameter es- 
timation an d retrieval methods (IMadhusudhan fe Seagerl 
' 2010all201lHMadhusudhan et al.ll2011aHLee et alJl2012t 
Line et all 12011a ). Potential drawbacks to these tech- 



niques are the necessary treatment of a three-dimensional 
(3D) system as a one-dimensional (ID) problem and the 
adoption of some built-in assumptions, such as the pres- 
ence and extinction characteristics of any clouds/hazes, 
or the selection of a limited number of molecules that 
are presumed to contribute to the observed behav- 
ior. Moreover, the best mathematical solution may 
not necessarily make physical sen se . As an exam - 
ple. IMadhusudhan fc Seagerl (|2009ft . iLee et all (|2012t ). 
and lLine et al.l (|2011blT all derive an unexpectedly large 
CO2 abundance for HD 189733b based on the dayside 
eclipse Hubb l e Spac e Telescope (HST) NICMOS data of 
iSwain et~a l. (2009aJ) that so far defies theoretical expla- 
nation, e ven when disequilibrium processes are consid- 
ered (see lMoses et al.l 120111: iLine et al.ll2010D . It is not 
clear at the moment whether the theoretical models are 
at fault (e.g., they could be missing some key mecha- 
nisms that enhance photochemical production of CO2) or 
whether there are some problems with the data and/or 
analysis (e.g., underestimated systematic error bars for 
the data, or missing opacity sources for the spectral mod- 
els). Although these techniques are powerful, it is always 
a good idea to follow up such studies with physically 
based forward models to confirm that the derived char- 
acteristics can be explained theoretically. 

One interesting and surprising result from these and 
other observational analyses is the suggestion of a low 
derived water abundance for several of the close-in tran- 
siting exoplanets examined to date, in comparison with 
CO and/or CH4, or in comparison with expectations 
based on solar-composition atmospheres in chemical 
equilibrium. The H 2 abundance h as been suggested to 
be ano malously low on WASP-12b ()Madhusudhan et"aL1 
I2011al see also iCowan et al.l I2012D , and this conclusion 
seems firm based on model comparisons with the Spitzer 



Space Telescope dayside eclipse data of iCampo et al.l 
(|2011[) . Water h as also been r e porte d to be low 
on HP 189733b dGrillmair et al.l [20071: ISwain et al l 
2009al: IMadhusudhan fc Seagerll2009l: iPesert et al.ll2009b 
Sing et all 120091 : [Line et al.l l20lTbl [Lee et al l 120121) 
on HP 209458b jSeager et al.l 120051: iRichardson et al l 
20071: ISwain et al.ll2009bt IMadhusudhan fc Seagerl 120091 



2010al) . and on XO-lb and CoRoT-2b flMadhusudhanl 
20121 ). The interpretation for these latter planets is not 



definitive given uncertainties in the thermal structure 
and other model parameters, and given that solutions 
with solar-like water abundances have been found that 
provide acceptable fits to the transit and eclips e data 
for these plan ets (e g. IBurrows et al.l 120051 l200l 1 20081: 
i Fortnev et "ail 120051 l2008al IFortnev fc Marievl 120071: 
Barman et~atl 120051: iBarmanl 120071 120081: iTinetti et al. l 
2007L 120101: iGrillmair et al.1 120081 Icharbonneau et al l 
20081: ISwain et al.1 l2008bt iMachalek et al.l 120081: 
Madhusudhan fc Seagerl 120091 l2010al iBeaulieu et al.l 
l2010t iKnutson et al.l I2012H 
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fit is often i mproved for t he case of 
abundances ( Madhusudhan 



because the 
low er water 
Seagen 120091 EfflTOal : 
IMadhusudhan et al.l l2011at IMadhusudhan! 12012ft . the 
suggestion of lower-than-expected H2O abundance 
should not be dismissed out of hand. 

Water is thermodynamically stable under a wide va- 
riety of temperature and pressure co nditions expected 
for hot- Jupiter atmospheres (e.g., IBurrows fc Shar 
1 19991: lLodders fe Feglevl 120021 : iHubenv fc Burrows! 1200 
iSharp fc Burrows! 12007). and due to the large cosmic 
abundance of oxygen, H 2 is expected to be one of 
the most abundant molecules on these transiting plan- 
ets if one assumes solar-like elemental abundances. Al- 
though H2O can be destroyed by ultraviolet photons, 
photochemical models predict that water will be effec- 
tively recycled in hot- Jupiter atmospheres (iMoses et al.l 
l20Tll ILine et~al|[20ll l20lTal) . keeping H 2 abundances 
close to thermochemical-equilibrium predictions. Extinc- 
tion from clouds and hazes could help mask the water 
absorption signatures, but in that case other molecu- 
lar bands would be affected as well, makin g the spec - 
trum appear more like a blackbody (e.g., iLioul l200"2h . 
If the interpretation is robust, one possible explana- 
tion could be that the atmospheres are relatively en- 
riched in carbon, with C/O ratios > 1, which leads 
to small H2O/ CO and H2O / CH4 ra tios in hot- Jupiter 
atmo spheres (ISeager et al.1 120051 : iKuchner fc SeageT 
20051: iFortnevI 120051: lLodders! 120101: IMadhusudhan et al 



2J)lla|bJ; [Madhusudhan 2012; Koppa rapu et al.ll201^ 



The influence of elemental abundances such as the 
bulk C/O and C /N ratios on exoplanet composition and 
spectra has not been as well studied as some of the 
other factors that can affect the observations — inves- 
tigators tend to simply assume solar-like elemental ra- 
tios in their models. This neglect is unfortunate, as 
giant planets can conceivably have elemental composi- 
tions different from their host stars, and the atmospheric 
C/O ratio can provide a particularly important clue 
for unraveling the formation and evo l utionary history of 
the p l anets (e.g.. iLunine et al.ll2004t lOwen fc Encrena: 
120061: Uohnson fc Estradal 120091; "iMousis et all 120111 
IMadhusudhan et al.ll2011bl : IOberg et al.l2011l ). The C/O 
ratio in a giant planet's atmosphere will reflect not only 
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the composition of the circumstellar disk in which the 
planet formed, but physical properties of the disk (such 
as the total disk mass and temperature structure), the 
planet's formation location within the disk (especially 
in relation to various condensation fronts), the rate of 
planet formation, and the evolutionary history of the 
planet and disk, including migration history and relative 
contribution of gas versus solid planetesimals in deliver- 
ing heavy elements to the atmosphere. 

Regardless of whether the planet formed and re- 
mained beyond what was once a putative "snow line" 
in the disk, like our own solar-system giant planets, 
or whether it migrated from behind the snow line 
or from elsewhere in the disk to its present short- 
period orbital positions, like the close-in transiting 
hot Jupiters, the C/O ratio can help unravel the 
conditions in the protoplanetary disk and the timing, 
history, and m echanism of th e planet's formation and 
evolu t ion (e.g.. | Larime nll975t IL arimer fc Bartholomayl 



197 ^ iLattimer et all I197S 



Stevenson fc Lunind 



Loddcrs & Fcglcv 
19981 ICvr etal 



" iLunine fc Stevenson! 
iPrinn fc FeglevI 
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Gautier et all l200ll iHersant et al 



2004 |2010T ILunine et all 12004 ; 



1993, 119951 119991: iNiemann et all 
1998, 1999: lOwen et all 11999: 
iKrot et afT [20001: iGaidosI 12 000; 

12004 120081; iLoddersI 
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Wong et al 

Cuzzi fc Zahnld 12004 iCuzzi et all 12005); lAlibert et al 

ICiesl 



Kuchner & Scager 



lOwen fc Encrena: 



iCiesla fc Cuzzi 



2006: iWoodenl 12008: D odsor 



Mousis et all I2009allbl. 



iGuillot fc Hueso 



.20111 

Madhusudhan et al.l 1201 lbl; IMadhusudhanl 
Naiita et aD 1201 It lOberg et al.ll20Tlt lEbel fc Alexander! 
20111: iKopparapu et al.ll2012h . 

The atmospheric C/O ratio itself can have a strong 
influence on the relative abundance of the spectro- 
scopically import ant H 2 Q, CH4, CO, CO 2 , C2H2, and 
HCN molecules (iLodders fc FeglevI [2001 iSeager et all 
20051: iKuchner fc Seagerl 120051: iZahnle et all 12009a!: 
Madhusudhan & S eagerll2010al: lL oddcrs 2 0101: ILine et all 
2010t IMadhusudhan et all l2011at IMadhusudhanl 120121; 
Kopparapu et al.l 120121 ). A photospheric value of C/O 
~ 0.55 seems to be the consensus for our Sun, despite 
some d isagreement over the exact absolut e C/H and O/H 
ratios (lAljende Prieto et all 120011 120021; lAsplund et all 
20051 l2009HGrevesseet all2007Ll2010HCaffau et alll2008l 
2010t lLodders et alll2009D . but possible departures from 
the canonical solar value or from the host-star's value 
need to be considered in atmospheric models. 

This need is spotlighted by the low derived H2O abun- 
dances for some exoplanets mentioned above, and in 
particular, by the recent i nference of a carbon - rich at - 
mospher e for WASP-1 2 b (IMadhusudhan et ail 1201 lal) 
In fact , IMadhusudhanl (|2012f ) and IMadhusudha n et all 
(|2011bf) suggest that super-solar C/O ratios could po- 
tentially explain not only the unusual eclipse pho- 
tometric band ratios of several transiting exoplanets 
(e.g., WASP-12b, WASP-14b, WASP-19b, WASP-33b, 
XO-lb, CoRoT-2b), but also the lack of an appar- 
ent thermal inversion in some highly irradiated hot 
Jupiters, due to the fact that TiO — one candidate ab- 
sorber that could cause a stratospheric thermal inversion 
(|Hubenv et alll2003tlFortnev et al.l[2006bl : lBurrows et all 



Robinson et al.l 



iBond et al.l 



2006) — does not form in significant quantities when at- 
mospheric C/O ratios > 1. To provide a good fit to 
the seco ndary-eclipse data, the carbo n-ric h models iden- 
tified by IMadhusudhan et al.l (|2011al ) and Madhus udhan! 
(2012) all require specific constraints on species abun- 
dances and thermal profiles, and it remains to be demon- 
strated that the derived abundances make sense from a 
disequilibrium-chemistry standpoint. 

We therefore develop equilibrium and disequilib- 
rium models to better quantify the effects of the 
C/O ratio on exoplanet atmospheric composition in 
a general sense, as well as to specifically check 
whether th e molecular abundances derived from analy - 
ses such as IMadhusudhan fc Seagerl (l2009ll2010 al| 201lL 
Madhusudhan et all (l2011afl, Madhus udhan! (2012). 
Lee et all <|2012f ) . and iLine et all (|2011bl ) are physically 



and chemically plausible for assumptions of either equi- 
librium or disequilibrium conditions. We will also check 
whether additional molecules that are typically not con- 
sidered in spectral models could be contributing to the 
observed behavior of transiting hot Jupiters. Our dise- 
quilibrium models include the effects of thermochemical 
and photochemical kinetics and vertical transport in con- 
trolling the vertical profiles for the major neutral carbon, 
nitrogen, and oxygen species through the entire relevant 
atmospheric column, from the equilibrium-dominated 
deep troposphere up through the quench regions for 
the major species, and on up t o the photochemistry- 
dominated upper strat o sphere (Visscher et alJ l2010bt 
iMoses et all 120101 l201lt iVisscher feMosesI 120111 ). Dif- 
ferences that occur as a result of the C/O ratio will be 
highlighted. 

2. MODELS 

We restrict our study to C-, N-, and O-bearing species 
(along with H and He) because these are the most cosmi- 
cally abundant of the reactive heavy elements, because 
they will be incorporated into the most abundant and 
most spectrally significant molecules in the infrared pho- 
tospheres of hot Jupiters, and because their kinetic be- 
havior is better understood than that of other elements 
like sulfur, phosphorus, alkalis, halides, or metals. The 
chemistry of other elements will be interesting in its own 
right, but given the lower cosmic abundances of these 
elements, their coupled chemistry with C-, N-, and O- 
bearing compounds is not expected to have a significant 
impact on the molecular abundances predicted here. One 
exception is the loss of ^20% of the oxygen if t he planet 
is coo l enough for silicate condensation to occur ([Loddersl 
2010), and we account for that loss in our model. 

We focus on transiting planets that are hot enough 
that CO rather than CH4 is expected to be the domi- 
nant carrier of atmospheric carbon because the compo- 
sition of these warmer planets is expect ed to be more 
sensi tive to changes in the C/O ratio (Madhusudhan 
I20"l2lh We ignore planets with well-defined strato- 
spheric ther mal inversions becau se we assume (perhaps 
naively; see iKnutson et al.l 120101 ) that Ti O is responsi- 
ble f o r upper-atmospheric heating (e.g.. | Hubeny et all 



ble to r upper- atmospheric heating (e.g.. Irlubeny et ai.l 
200 1 iFortnev et al.l l2006bl l2008at iBurrows et all I2006L 



2007 , 120081) . and thus these planets may be expected t o 



have solar-like C/O ratios (|Madhusudhan et alll2011bl) . 
The disequilibrium c hemistry of such p l anets has already 
been explored (e.g., IMoses et all 1201 it ILine et all 120101 : 
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IZahnle et al.ll2009bt iLiang et alJl200ll200l . 

For our equilibrium calculati ons, we use the 
NASA CEA code developed by iGordon fc McBridd 
(Il994fl. with thermodynam ic parameters adapte d from 
Guryich et all (11989-1994D the JANAF tables (|Chasd 
19981 ). IBurcat fc Ruscid (|2005f ) . and other literature 



sources. For the disequilibrium calculations, we use the 
Calte ch/JPL ID photo chemistry /diffusion code KINET- 
ICS (|Allen et al.l Il981f) to solve the coupled continuity 
equations that control the vertical distributions of tropo- 
spheric and stratosp heric species in th e case of incident 
e .g.. lYung et al.l I1984J: IGladstone et al 



UV irradiation 



i .e., I Yung 

19961 : iMoses et al.ll2000allbl. 2005. 201 0ll2011l;lLiang et all 
20031 l2004tlLine et al.ll20ldl2011at" Tvisscher et al I l2010bt 
Visscher fc Mosesll201lD . " 



2.1. Photochemical and Thermochemical Kinetics and 
Transport Model 

When transport time scales are shorter than chem- 
ical kinetics time scales, the atmospheric composi- 
tion can depart from strict thermodynamic equilib- 
rium (e.g., Prinn fc Barshavl[T977tlLewis fc Fegle^HHl 
IVisscher fc Mosesl 1201 If) , making thermochemical equi- 
librium assumptions no longer valid for spectroscopic 
modeling. This so-called "transport-induced quench- 
ing" disequilibrium mechanism will affect the observed 
composition. The consequences f or extrasolar gi- 
ant p lan ets have been explo re d by ILodders fc Feglevj 



(2002), 



Cooper fc Showman] (|2006f) . Visscher et al 



Fortnev et al.l 



2007 V ILine et al l (1201 



(2006a), 
I2011a|) 



Hubenv fc Burrows! 



IMoses et al.l (120111). 



Vis scher fc MosesT 



, and IMadhusudhan fc Seagen 
(|2011l) . The atmospheric composition can also be 
altered by the intense ultraviolet radiation from the 
Absorption of UV photons by at- 
species can lead to dissociation of the 
and the resulting species can react to 
disequilibrium constituents. Models of 
of transitin g hot Jupite rs have 

aSS). 

Id] 



host star, 
mospheric 
molecules, 
form new 
the photochemistry 
been p resented bv ILiang et all (120031 1200 
(12004ft. [Garcia Munod (l2007D. IKoskinen et 



Zahn leetall (l2009allbT) . ILine et al.l (120101. 



Mos es et al.f(l2011l). iMiller-Ricci Kempton et all (201 




Kopparapuet al.l (|2012h . and iVenot et al.l (1201a) . The 
bottom line from the most recent models is that pho- 
tochemical processes will affect the composition of the 
middle atmospheres (i.e., the infrared "photospheres") 
of all but the hottest of the hot Jupiters, leading to the 
removal of some expected species like CH4 and NH3 from 
high altitudes and to the production of some other inter- 
esting spectrally active constituents like atomic C, O, N, 
and molecular HCN, C2H2, and in some instances, CO2. 

For our disequilibrium models, vertical transport is as- 
sumed to occur via molecular and eddy diffusion. The 
eddy diffusion coefficient K zz , which is a free parame- 
ter in our models, provides a measure of the strength of 
atmospheric mixing via convection, atmospheric waves, 
and eddies of all scales. Although the vertical winds de- 
rived fr om general-circu l ation models (GCMs) such as 
those of iShowman et all (I2009D can be used as a guide 
for defining K„ (see IMoses et al.l 120 111) , we simply as- 
sume an altitude-independent (constant) value for K zz 
for the models presented here. We allow the solutions for 
the vertical species profiles to reach a steady state. The 
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Fig. 1. — Temperature profiles (solid lines, as labeled) adopted 
for our dayside atmospheres, s elected from models present ed in 
Madhusudhan & Scagcr (2009), Madhusud han et alj j |2011a! ), and 
Madhusudhan ( 2 0121) . The gray dot-dashed lines illustrate the con- 
densation curves for enstatite (MgSiOa) and forsterite (Mg2Si04), 
and the gray dashed lines indicate equal-abundance curves for 
N2-NH3 and CO-CH4, all for assumed solar-composition, solar- 
metallicity atmospheres. The more reduced components (CH4 and 
NH3) dominate to the left of the gray-dashed curves, and the more 
oxidized species (CO and N2) dominate to the right. Note that 
because our solutions do not necessarily favor solar-composition 
atmospheres, the gray curves will need to be shifted around in re- 
sponse to specific model assumptions; however, CO and N2 are 
expected to be the dominant carbon and nitrogen components in 
the photospheres (a few bar to ~0.1 mbar) of all these planets for 
all our assumed bulk elemental compositions and metallicities. A 
color version of this figure is available in the online journal. 

models all have 198 vertical grid points spaced equally 
in log(pressure), typically ranging from the deep tropo- 
sphere (^1 kbar) to well into the thermosphere (pressures 
<10 -10 bar). The lower boundary is chosen at a hot- 
enough temperature (typically > 2400 K) to encompass 
the quench levels for the various major disequilibrium 
species (e.g.. IMoses et al.ll20H"HVisscher fc Mosesll2011l) . 
and the upper boundary is chosen at a high-enough al- 
titude such that the atmosphere is optically thin to the 
UV radiation that dissociates the key photochemically 
active species in the model. Note, however, that we do 
not consider ion chemistry or other processes that are 
important in the thermosphere and upper stratosphere, 
and our results at pressures less than ~ 1 /ibar should not 
be taken seriously. For discu ssions of the r mospheric com- 
positio n and chemistry, see [Yelld (120041) , iGarci'a Munod 
(|2007t ). and IKoskinen et al.l (|2010| ). Plane-parallel ge- 
ometry is assumed for the radiative calculations in the 
photochemical model, and multiple Rayleigh scattering 
by H2 and He is considered using a Feaut rier radiative- 
transfer method (Mi chelangeli et alJll992f ). 

The thermal structure is an input to the KINETICS 
code; that is, temperatures are not calculated self consis- 
tently. We use specific profiles from the se conda r y-eclips e 
analyses of IMadhusudhan fc Seagerl (120091 [2 010a). 
IMadhusudhan et al.l (|2011afl . and iMadhusudhanl ( 2012) 
to define the thermal structure (see Fig. [IJ. The 
dayside hemisphere of the planet contributes to the 
secondary-eclipse observations, and in fact the hottest 
regions on the dayside will disproportionately domi- 
nate the planetary flux; therefore, our adopted pro- 
files are typically hotter than those derived from ID 
radiative-equilibrium models with efficient heat redis- 
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tribution (e. g.. iBurrows et all 12006. 20081: iFortnev et all 
I2006U [2010HBarmanll2008D . In principle, it would seem 
better to use results from the 3D GCMs to define our 
thermal structure; however, 3D models are not available 
for most of the planets we are investigating, and the com- 
plexity of the real atmospheres will make it difficult to 
realistically simulate all the chemical and physical mech- 
anisms that affect the spectrally active molecules and 
other opacity sources in the GCMs. We therefore rely 
on the statistical models listed above, which at least sat- 
isfy energy balance and provide an adequate reproduc- 
tion of the observed eclipse data, to define our thermal 
profiles. Since the quenching of the nitrogen species oc- 
curs at high temperatures and pressures deeper in the 
planet's troposp here than is accounted fo r with t he ther - 
mal profiles of IMadhusudhan fc Seaeerl (120091 l2010aft . 
IMadhusudhan et all (|2011al ). and iMadhusudhanl (|2012( i. 
we extend these profiles to deeper pressures assuming a 
somewhat arb i trary a diabat based on the ID models of 
IFortnev et al.l (|2008al) . At the low-pressure end, we ex- 
tend the profiles upward nearly isothermally (see Fig.[T]). 
The real upper atmospheres wil l transition into a high- 
temperatu r e therm osphere, but iMoses et al.1 (|2011| ) and 
I Line et all (|201 laf ) demonstrate that the adopted ther- 
mospheric structure has little effect on the results for the 
observable portions of the troposphere and stratosphere. 
Assuming these fixed temperature-pressure profiles, we 
then use the hydrostatic equilibrium equation to define 
the altitude grid. 

The reaction list is taken from lMoses et all (|2011f ). and 
the reader can refer to that paper for further details. The 
reactions are fully reversed assuming thermodynamic re- 
versibility, which allows the model atmosphere to reach 
thcrmochemical equilibrium when kinetic reaction time 
scales are much shorter than transport time scales. Our 
nominal 1 x solar composition atmospher es assume the 
protos olar composition from Table 10 of lLodders et al.l 
(2009); however, because our model does not include the 
rock-forming elements, we assume that 21% of the oxy- 
gen is tied up with condensed silicates and metals and 
is therefore not available for gas-phase chemistry (e.g., 
IVisscher fc Feglevl [200l ILoddersI l20Tol) . Our C/O ra- 
tios therefore refer strictly to the gas-phase C/O ratio 
after any silicate and metal-oxide condensation has oc- 
curred and not to the bulk planetary C/O ratio. Note 
that when we alter the C/O ratio above solar, the oxygen 
abundance is kept fixed while the carbon abundance is 
adjusted. Thermochemical equilibrium abundances from 
the CEA code are used to define the initial species pro- 
files, and zero flux boundary conditions are adopted at 
the top and bottom boundaries. Again, as discussed in 
IMoses et al.l ()2011[ ). we expect the atmospheres of these 
planets to be actively escaping, but our upper boundary 
condition has no effect on the stratospheric and tropo- 
spheric profiles unless the hydrodynamic wind extends 
down below the thermosphere, at which point our hy- 
drostatic models will be no longer valid. 

For the incident stellar ultraviolet flux that drives the 
photochemistry, we use HST/STIS spectra from epsilon 
Eri (a K2 V star) a nd Chi 1 Ori (a GO V star) from the 
CoolCAT database (|Aversll2005t) for the 1150-2830 A re- 
gion to mimic spectra from HD 189733 and WASP-12, 
respectively, after adjusting for the appropriate distance 



scalings. For wavelengths above and below this range, 
we scale the solar spectrum to match the expected mag- 
nitude of the flux for the appropriate stellar type and 
orbital distance. For XO-lb and CoRoT-2b, we sim- 
ply adopt the solar spectrum, adjusted to the appropri- 
ate orbital distance. The stellar zenith angle is fixed at 
48°, which is appro priate for secondary-eclipse conditions 
(jMoses et alJl20Tl . 

2.2. Spectral Models 

Given the atmospheric thermal and chemical profiles, 
as described above for each planet under consideration 
(namely XO-lb, HD 189733b, CoRoT-2b, and WASP- 
12b), we use a radiative transfer code to generate a 
corresponding model spectrum for comparison with ob- 
servations. We compute emergent thermal spectra us- 
ing the exoplanetary atmosph eric modeling program of 
IMadhusudhan fe Seagerl (|2009f ). The code computes ID 
line-by-line radiative transfer in a plane-parallel atmo- 
sphere under the assumptions of local thermodynamic 
equilibrium (LTE) and hydrostatic equilibrium. Line-by- 
line molecular absorption due to H2O, CO, CH4, CO2, 
NH 3 , C2H2, HCN, and H 2 -H 2 collision-induced absorp- 
tion (CIA) is considered. The molecular linelists for H2O, 
CO, C H4, and NH3 are obtained from iFreedman et al.l 
(2008) and references th erein, and the HITRAN database 
(jRothman et al.1 120051 12009( 1 . T he linelists for CQ 2 , 
C2H2, and HCN are obtained from iWattson fc Rothmanl 
(p8l . lRothman et al.l (|200l . and Harris et al.l (|2008). 
respectively. The CIA opacities are o btained from 
iBorvsow etaLl (|1997h and iBorvsowl (|2002( 1. The emer- 
gent spectrum is calculated assuming the planet is at 
secondary eclipse, i.e. that the integrated flux from the 
entire day side atmosphere of the is observed. In order 
to co mpute planet-star flux rat ios, a Kurucz model spec- 
trum (C astelli fc Kuruc3l2004fl is adopted for the given 
stellar parameters. While comparing our models with ob- 
servations of a given planet, we integrate the model spec- 
tra over the instrumental bandpasses for which data are 
avail able for that system and compute a goodness of fit 
(e.g.. IMadhusudhan fc SeageHl20"o1l IMadhusudhan et al.l 
I2011al) . 

3. RESULTS 

Assuming a "nominal" temperature-pressure profile se- 
lected for each of the exoplanets under consideration (see 
Fig. [T]), we investigate how the atmospheric chemistry 
changes for different assumptions of the atmospheric 
C/O ratio and metallicity. We first calculate thermo- 
chemical equilibrium constituent abundances as a func- 
tion of pressure for C/O ratios ranging from 0.1 to 1.9, 
for a variety of assumed bulk atmospheric metallicities. 
Then, we select one or more representative models with 
different metallicities and C/O ratios to investigate pos- 
sible disequilibrium chemistry in more detail with our 
thcrmochemical and photochemical kinetics and trans- 
port models. Disequilibrium-model solutions will be 
presented for carbon-to-oxygen ratios of both solar-like 
(C/O ~ 0.5) and carbon-rich (C/O > 1) atmospheres. 
We then create synthetic spectra from the model results 
and compare these with the secondary-eclipse data to 
determine if we can find disequilibrium model solutions 
that provide an acceptable fit to the data. Our result- 
ing models will be non-unique. Note that because haze 
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extinction is more likely to affect the transit observa- 
tions (due to the long slant path being able to amplify 
extinction from layers that would otherwise be optically 
thin in the vertical), we restrict our data comparisons to 
the secondary-eclipse observations. In addition, because 
Spitzer photometric data are available for all the planets 
in question, we focus on the Spitzer dataset. 

For any given metallicity, thermal structure, and fi- 
nite list of molecules that contribute to the spectral sig- 
natures, the abundances of the major gas-phase opacity 
sources required to fit the secondary-eclipse data are well 
constrained, and there will be a correspondingly narrow 
range of C/O ratios (if any) that will produce these de- 
sired abundances. However, the best-fit C/O ratio will 
change for different assumptions about the thermal struc- 
ture and metallicity. Although we examine a few differ- 
ent scenarios in this study, we do not explore the full 
range of possible parameter space. Our main goal is to 
illustrate how the C/O ratio affects the composition of 
highly irradiated hot Jupiters. Our purpose in devel- 
oping individual models with acceptable spectral fits is 
to determine whether atmospheres with chemical abun- 
dances that are theoretically consistent with "reason- 
able" metallicities and thermal structures can be derived 
that reproduce the observed secondary-eclipse data, and 
in particular, whether models with high C/O ratios are 
reasonable for the exoplanets that have been reported to 
have relatively low H 2 abundances. Although we fo- 
cus on four specific transiting planets in this paper (XO- 
lb, HD 189733b, CoRoT-2b, WASP-12b), the planets we 
have selected span a wide ran ge of apparent dayside tem- 
peratures from -1300-2900 (jCowan fc Agol l20lTl ) and 
can thus serve as proxies for a variety of hot Jupiters. 

3.1. HD 189733b Results 

The moder ate-temperatur e exop lanet HD 189733b, 
discovered by iBouchv et al.l (|2005|), resides in the up- 
per part of the 01 regime of the [Madhusud hanl (|2012f ) 
classification scheme, in which the C/O ratio is less 
than unity (as is favored from s tatistical analyses of the 
HD 189733b observational data; Madhus udhan fc Seagerl 
2009), and temperatures are warm enough that CO dom- 
inates over CH4 in the observable portion of the atmo- 
sphere (see Fig. [TJ, yet temperatures are cool enough 
that Ti will be tied up in condensates rather than 
as gas-phase TiO. Various analyses suggest that HD 
189733b does not have a significant stratospheric ther- 
mal inversion (jFortnev et~ . 2008a; Burrows et al.ll2008t 
Barman 2008; Charbonneau et al. 2008; Grillmair et al. 
2008t iSwain et al.l l2009at IMadhusudhan fc Seagerl 120091: 
Showman et al.ll2009fl . ~ 

3.1.1. HD 189733b Chemistry 

Figure [2] illustrates how the thermochemical equilibri- 
um composition of HD 189733b varies as a function of 
the atmospheric C/O ratio for our assumed nominal day- 
side thermal structure and an assumed lx solar metal- 
licity. For a solar-like C/O ratio of —0.5, CO and H2O 
are the most abundant heavy molecules, CO2 is much 
less abundant at the < 1 ppm level, CH4 and NH3 are 
well below ppm levels in thermochemical equilibrium at 
pressures less than 1 bar, and HCN and C2H2 are trace 
constituents in equilibrium such that they would have 



no influence on the eclipse spectra. Although transport- 
induced quenching could lead to CH4 (and NH 3 ) mole 
frac tions greater than 1 ppb at pressures less t han 1 bar 
(see lMoses et aUteOlUlVisscher fc Mosesll2011h . the CH 4 
mole fraction is never likely to rival that of CO for solar- 
like C/O ratios, unless effective eddy diffusion coefficients 
are gre ater than 10 11 cm 2 s" 1 in the CO-CH4 q uench 
region (|Moses et al.l 12011 fvlsscher fc Mosedl20lTI ). The 
relatively low derived CH4 abundances from the tran- 
sit and eclipse observations (ISwain et al.l l2008bl. [2 009a; 
Madhusudhan fc Scagcr 2009t lLee et al.ll2012HLine et al.l 
201 lbl ) suggest that atmos pheric mixing is no t that vig- 
orous on HD 189733b (see Ivisscher fc Mosedl20H . and 
CO mole fractions are likely greater than those of CH4. 

On the other hand, when the C/O ratio is greater 
than unity, Fig. [5] shows that the equilibrium methane 
abundance can approach or even exceed that of CO; 
other oxygen-bearing species like H2O and CO2 become 
much less abundant, and previously unimportant carbon- 
bearing species like HCN and C2H2 become abundant 
enough to influence spectra. At dayside photospheric 
temperatures and pressures (—1800-900 K, — l-to-10 -4 
bar), CO is a major constituent for all assumed C/O 
ratios. When C w O, available carbon and oxygen will 
bond together to predominantly form CO; when excess 
O is available (i.e., C < O), the extra O goes into the 
formation of H2O, and to a lesser extent CO2; when ex- 
cess C is available (i.e., C > O), the extra C goes into the 
formation of CH4, and to a lesser extent HCN and C2H2. 
When the atmospheric C/O ratio is very near unity, small 
changes in the C/O ratio lead to large changes in the 
species abundances. 

Wor king solely from the Madhusudha n fc Seagerl 
(2009) statistical constraints on the ensemble of mod- 
els that provide the best fit to the HD 1 8 9733b Spitzer 
seco ndary-eclipse data of iKnutson et al.l (|2007l [2009b) 
and iCharbonneau et all (|2008h (see the red horizontal 
segments in Fig. [2]), it is clear that equilibrium models 
with C/O ratios < 1 will be required to provide a de- 
cent fit to the o bservations. Here we u s e the £ 2 < 1 
constraints from IMadhusudhan fc Seagerl <|2009f ) . where 
i 2 = X 2 /N bs, with N b s — 6 being the number of pho- 
tometric data points available for the goodness-of-fit cal- 
culations. The CH4 and CO2 constraints seem to favor 
low C/O ratios (perhaps even C/O < 0.1), whereas the 
H2O constraints can be satisfied with higher C/O ratios. 
In reality, the water abundance has the biggest influence 
on the spectral behavior. For our nominal temperature 
profile and an assumed lx solar metallicity, H2O mole 
fractions that fall near — lxlO -4 provide a significantly 
better fit to the data than mole fractions even a factor 
of a few away from the middle of the H2Q-constra int 
range provided by the Madhusudha n fc Seagerl (|2009| ) £ 2 
< 1 criterion. It turns out that for a given temperature 
profile and metallicity, there is a very narrow range of 
C/O ratios that provides the necessary best-fit H2O mole 
fraction. For our assumed conditions, a C/O ratio near 
0.9 provides the desired equilibrium water abundance 
of -lxl0~ 4 . On the other hand, a C/O ratio of -0.9 
also implies a CH4 abu ndance at the extreme h i gh en d 
of what is allowable by IMadhusudhan fc Seagerl (2009) , 
and implies a very low CO2 abundance that appears to 
fall well outside the £ 2 < 1 range for CO2. Keep in 
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Fig. 2. — Mixing-ratio profiles for H2O, CH4, CO, CO2, HCN, and C2H2 on HD f89733b from our thermochcmical-equilibrium models 
with assumed lx solar metallicity but different assumed C/O ratios (ranging from 0.1 to 1.9, incrementing by 0.1 — dashed black lines 
for C/O < 1, solid lines for C/O > 1). The colored lines represent disequilibrium-chemistry results for a model with lx solar metallicity 
and a C/O ratio of 0.5 (purple), a model with lx solar metallicity and a C/O ratio of 0.88 (green), and a model with 3x solar metallicity 
and a C/O ratio of 0.96 (blue). All disequilibrium models have an assumed moderately- low eddy diffusion coefficient of K zz = 10 8 cm 2 
s . The horizontal line segments (red) represent the abundance constraints d erived by Madhusudhan & Seager (2009) from Spitzer IRAC 
and MIPS observations IIKnutson et al.il2007U2009bl ; ICharbonneau et al.ll200l ). Note that IMadhusudhan He Seagerl H2009I 1 did not consider 
HCN and C2H2 in their analysis, so abundance constraints for these molecules are not available. A color version of this figure is presented 
in the online journal. 
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Fig. 3. — Mole-fraction profiles for several important species 
(as labeled) in a thermochemical and photochemical kinetics and 
transport model for HD f 89733b. Model assumptions include lx 
solar metallicity, a C/O ratio of 0.88, and K zz = 10 s cm 2 s -1 
(see also the green curves in Fig. [2J. The dotted lines show the 
thermochemical-cquilibrium abundances for the species, using the 
same color coding. A color version of this figure is available in the 
online journal. 

mind, however, that the CO2 ab undance is essentially un- 
const rained at the £ 2 < 2 level (Madhus udhan fe Seagerl 
2009), and the degeneracy between the contributions of 
CO and CO2 in the IRAC channels could be compli- 
cating the identification and quantification of the CO2 
abundance. 

Can disequilibrium chemistry help improve the com- 
parison with observational constraints? Figure [2] also 
shows results from a disequilibrium model with an as- 
sumed lx solar metallicity and atmospheric C/O ratio 
of 0.88 (the green curve). For this model, also shown 
in Fig. |SJ the H2O abundance essentially follows the 
chemical equilibrium curve until it reaches high-enough 
altitudes that molecular diffusion begins to control the 
profile, at altitudes well above the infrared photosphere 
(located at a few bars to ~0.1 mbar). Photochemistry 
and transport-induced quenching therefore do not change 
the favored C/O ratio of ~0.9 for a lx-solar-metallicity 
atmosphere and our nominal temperature profile. We 
also find that photochemical production of CO2 from CO 
and H2O chemistry results in little net additional CO2, 
except at very high altitudes, beca use the CO2 is abl e 
to recycle back to CO and H2O (see M oses et alj|20lil ). 
In this regard, the results from our fully-r eversed reac- 
tion m echanism are consistent with those of iMoses et al.1 
(I20T1 . who find that C0 2 remains in equilibrium with 
H 2 and CO throughout much of the day side atmo- 
spheri c column, in contrast to the results of iLine et al.l 
(2010), in which an incompletely reversed reaction mech- 
anism allows C0 2 abundances to be enhanced by a fac- 
tor of ~2 due to photochemistry. Even with a factor-of- 
two photochemical enhancement, the CO2 mole fraction 
would not be large e nough in the C/O = 0.88 case to 
fall within the desired IMadhusudhan fe Seagerl (|2009t l £ 2 
< 1 range for CO2 . For the case of methane, transport- 
induced quenching slightly enhances the mole fraction as 
compared with equilibrium in the ^l-to-10~ 5 bar region, 
and only above the infrared photosphere do disequilib- 
rium processes reduce the CH4 abundance; thus, disequi- 
librium chemistry does not improve the CH4 abundance 



compa risons with regard to the Ma dhusudhan fe Seagerl 
(2009) £ 2 < 1 range for methane. Given that water dom- 
inates the spectral behavior in the mid-infrared, it may 
not matte r that the other species do no t fall precisely 
within the Madhusudhan & Scager (2009) desired ranges 
(see Section l"3. 1 . 21 below) . 

The disequilibrium chemistry of HD 189733b with a 
so lar-like atmospheric C/O ratio is discussed extensively 
in IMoses et al.l (|2011f) , and we will not repeat that dis- 
cussion here (see also ILine et aLl 12010( 1 . The main re- 
action mechanisms controlling the composition are very 
similar in an atmosphere with C/O = 0.88 versus 0.5, 
but the resulting abundances can vary significantly be- 
tween the two cases because of the different equilib- 
rium stabilities. For instance, the water abundance is 
higher in the C/O = 0.5 disequilibrium model than in the 
C/O = 0.88 model purely because there is more oxygen 
available once carbon sequesters an amount of oxygen 
need to form CO, and because photochemical processes 
do not efficiently remove H2O and convert it to other 
oxygen-bearing species. The dominant H2O recycling 
mechanism o perat ing in our models is scheme (10) in 
IMoses et al.l (|2011l ). Because our no minal thermal struc- 
ture is warmer than that assumed in IMoses et al.l (|2011[ ) 
(and our adopted K zz is smaller), our que nched abun- 
dance s of CH4 and NH3 are smaller than in IMoses et al.l 
(|2011f) . Note from Fig. [2] that CH 4 quenches at about 
the same pressure level for both the C/O = 0.88 model 
(green) and C/O = 0.5 model (purple). The quench 
point is under th e control of reaction scheme (2) of 
IMoses et al.l (|20T1 . with OH + CH 3 -> CH 2 OH + H 
being the rate-limiting step for the quenching of CH4 — > 
CO conversion under the ther mal-structure condition s 
of our nominal model (see also IVisscher fc Mosesl[20iTI) . 
Because the equilibrium CH4 abundance is larger at that 
quench point in the C/O = 0.88 model, the disequilib- 
rium quenched methane abundance is larger with C/O = 
0.88 than with C/O = 0.5. The larger quenched methane 
mole fraction then enhances the effectiveness of mecha- 
nisms that con vert CH4 to hydro carbons like C2H2 (see 
scheme (12) of Moses et al. 2011), and the total column 
abundance of C 2 H 2 through the middle atmosphere is 
significantly greater for the C/O = 0.88 model. 

A similar result occurs for the HCN abundance at low 
altitudes, but not at high altitudes. Hydrogen cyanide 
can continue to remain in equilibrium with CO, CH4, and 
NH3 in the lower and middle stratosphere after these 
molecules quench. The abundance of molecules that 
do not contain carbon and oxygen, like NH3 and N2, 
tend to be relatively unaffected by the C/O ratio. How- 
ever, because the CO and quenched CH4 abundances are 
larger for greater C/O ratios, processes that kinetically 
convert CO and NH3 into HCN (see scheme (8) from 
IMoses et all 120111 ) or that conve rt CH 4 and NH 3 into 
HCN (see schemes (7) and (14) in lMoses et"aT1l201lH are 
more effective with high C/O ratios than low C/O ra- 
tios. That leads to a much larger HCN abundance at 
low- and mid-stratospheric altitudes above the quench 
point with the C/O = 0.88 model as compared with the 
C/O = 0.5 model. However, the upper atmospheric peak 
abundance of HCN is similar in both models due to the 
fact that CO and N2 photochemistry are responsible for 
the high-altitud e production (e.g. , see the first scheme 
in section 3.5 of IMoses et~al"1 120 111) . High-altitude HCN 
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Fig. 4. — Secondary-eclipse spectra for HD 189733b (in terms of 
the ratio of the flux of the planet to that of the host star). The solid 
lines represent synthetic spectra from the results of our disequilib- 
rium chemistry models with lx solar metallicity and K zz = 10* 
cm 2 s — 1 , for C/O ratios of either 0.88 (green) or 0.5 (purple). The 
blue circle s with error bar s repre s ent the Spitz er IRAC and MIPS 
data from [Knutson ct al. (2007, 20092) and ICharbonnca u et all 
(2008). The small bl ack circles with err or bars represent the 
HST/NICMOS data of ISwain et"all H2009af ). and the gray circles 
with e rror bars represent the Spitzer IRS data of Grillmair et al. 
(2008). The circles without error bars represent our C/O = 0.88 
(green) or C/O = 0.5 (purple) model results convolved over the 
Spitzer IRAC and MIPS instrument response functions. Note the 
increased absorption for the purple C/O = 0.5 model in all the 
Spitzer channels; the C/O = 0.88 model fits the MIPS and IRAC 
data better than the C/O = 0.5 model for our particular assumed 
metallicity and thermal structure. The dashed lines are blackbody 
curves for assumed temperatures of 900 K (bottom), 1400 K (mid- 
dle), and 1800 K (top). A color version of this figure is available 
in the online journal. 



production occurs through the reaction N + OH — >• NO 
+ H, followed by NO + C -)• CN + O, followed by CN 
+ H 2 -> HCN + H. The N and C required for this mech- 
anism are derived from N 2 and CO photolysis. Both 
these molecules become optically thick to the photolyz- 
ing UV radiation (and therefore self shield) at high al- 
titudes, so this mechanism only works within a limited 
altitude region, but the HCN production rates are sub- 
stantial. Because the CO and N 2 column abundances at 
optical depths less than a few are the key defining values 
for this process to work (which does not change between 
models), the total production rates from this process are 
similar in our C/O = 0.5 and C/O = 0.88 models. Hy- 
drogen cyanide is therefore a key photochemical product 
regardless of the C/O ratio, although the total column 
above ~1 bar will be larger for higher C/O ratios due to 
the thermoc hemical kinet i cs sch emes (7), (8), and (14) 
mentioned in Moses et al. (20 111) ■ 

In both models, CO and C0 2 track equilibrium profiles 
because kinetic recycling is efficient, and there are fewer 
"leaks" into or out of the cycles compared with CH 4 , 
NH 3 , HCN, an d C 2 H 7: hydrocarbo ns. The CO recycling 
is discussed in iMoses et al.l (|2011l ) . Carbon dioxide is 
recycled through the forward-reverse reaction pair OH 
+ CO ^> C0 2 + H, which is effective throughout the 
stratosphere for both models. Some net photochemical 
production for C0 2 does occur at high altitudes, but 
it does not change the total column abundance in the 
photosphere. 



3.1.2. HD 189733b Spectra 

Figure [4] shows synthetic spectra from our disequilib- 
rium models, in comparison with the Spitzer IRAC and 
MIPS secondary-eclipse data for HD 189733b. The plot 
shows spectra from disequilibrium models with assumed 
C/O ratios of 0.88 and 0.5 (the green and purple curves 
from Fig. [2]). These two chemical models have identical 
assumptions other than the C/O ratio. As is obvious 
from the plots, the model with C/O = 0.88 provides a 
statistically better fit to the data than the more solar-like 
C/O = 0.5 model for our assumed 1 x -solar- metallicity 
atmosphere with our adopted nominal thermal profile 
(i.e., £ 2 < 1.9 for the C/O = 0.88 model versus £ 2 < 
5.0 for the C/O = 0.5 model). It is also apparent from 
the plot that the fit could have been improved with some- 
what less methane (note the excess model absorption at 
3.6 and 8.0 microns, where CH4 contributes opacity) and 
perhaps more C0 2 (note the fit at 4.5 microns, where 
CO and CQ 2 contribute opacity), as was suggested by 
the lMadhusudhan fe Seagerl (|2009f ) £ 2 < 1 ranges. How- 
ever, the fit from the disequilibrium chemical model with 
C/O = 0.88 is within 1.4-sigma in all the Spitzer bands, 
on average, implying that theoretically plausible models 
with carbon-enhanced atmospheres (compared with the 
solar C/O ratio) can be consistent with the HD 189733b 
Spitzer secondary-eclipse observations. 

It must be kept in mind that the fit is not unique. 
More solar-like C/O ratios may also be consistent with 
the data, but in that case, Fig. [4] demonstrates that 
the atmosphere must be much warmer than we are as- 
suming in our nominal model, as our current combina- 
tion of thermal structure and its resulting chemically- 
constrained water abundance in the C/O = 0.5 model 
leads to too much absorption in many of the Spitzer 
bandpasses. However, our temperatures are already at 
the maximum of what has been predicted from self- 
consistent I D models with no energy redistribut ion from 
the dayside (lBarmanll2008l iBurrows et al.ll2008t see also 
iFortnev et al.ll2006bf ) or from the dayside-ay erage atmo- 
sphere in 3D models ()Showman et al.l I2009T ) ■ so there is 
a limit to how low the C/O ratio can go to fit the data 
while still maintaining a plausible thermal structure. 

The solution also depends on assumptions about atmo- 
spheric metallicity. For a fixed thermal structure, differ- 
ent metallicities require a different C/O ratio to maintain 
the desired H 2 mole fraction to compare well with ob- 
servations (see Fig. [5]). For example, if the atmospheric 
metallicity were 3x solar, the C/O ratio would need to 
be ~0.96 for the H 2 mole fraction to be near 10 -4 in 
equilibrium in the photospheric region of HD 189733b 
for our nominal thermal structure (see also the disequilib- 
rium model represented by the blue curve in Fig. [2]). The 
greater the metallicity, the closer to unity the C/O ratio 
must be to remain consis t ent wi th the spectral analysis of 
Madhusu dhan fe Seagerl (|2009D (i.e., super-solar metal- 
licities tend to require super-solar C/O ratios to remain 
consistent with observations, for atmospheric tempera- 
tures at or below those represented by our nominal ther- 
mal profile). For subsolar metallicities, more solar-like 
C/O ratios could pr oduce the desired H 2 Q abundance. 

As was shown by iMadhusudhanl (|2012l ) , the relative 
abundance of methane and water is a sensitive indicator 
of the C/O ratio in hydrogen-dominated atmospheres, 
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Fig. 5. — The mole fraction of H2O (Left), CO2 (Middle), and the CH4/H2O ratio (Right) as a function of the C/O ratio and metallicity, 
assuming thermochemical-equilibrium conditions at 1345 K and 80 mbar (i.e., within the photosphere of our nominal model for HD 
189733b). The shaded regions represent the £ 2 < 1 observational constraints from Madhusudhan & Seaecr (2009). The CH4/H2O ratio is 
very sensitive to the C/O ratio, and the CO2 and H2O mole fractions are very sensitive to both the atmospheric metallicity and C/O ratio. 
The planetary spectrum itself is very sensitive to the water abundance, and the red curve in the figure on the left shows our best-fit H2O 
model fraction for our nominal thermal structure. The observations indicate C/O < 1 for HD 189733b, but not all of the observationally 
allowable phase space is chemically plausible. The higher the atmospheric metallicity, the closer to unity the C/O ratio must be to provide 
the best fit to the Spitzer secondary-eclipse spectra. A color version of this figure is available in the online journal. 



particularly when photospheric temperatures are warm 
enough that CO dominates over methane as the ma- 
jor carbon carrier. A CH4/H2O ratio less than unity 
indicates a C/O ratio l ess th an 1 (see Fig. [5]). The 
iMadhusudhan fc Scagei] (|2009T ) analysis firmly suggests 
a C/O ratio < 1 for HD 189733b; however, just how 
far below unity the C/O ratio can be to remain consis- 
tent with observations depends on metallicity. Figure [5] 
demonstrates t hat while CH4/H2O ratio s < 10 ~ 5 are al- 
lowable by the IMadhusudhan fe Seagerl (|2009l ) spectral 
analysis, such values are not physically plausible unless 
the C/O ratio is strongly subsolar (i.e., < 0.1) and/or 
metallicities are greatly supersolar. 

If the CO2 abundance were well constrained from 
observations, one could also constrain the overall 
atmospheric metallicity due to the well-known sen- 
sitivity of CO2 to me tallicity (Fig. [5l see also 
Fodders fc Feglevl 120021 iFortney et all |200" 
2008U IBurrows et all 120061 120081 I Zahnle et all 



2006b 



2009a, 



Line et al.ll2010L IMoses et al J 1 20 111 ). However, using the 



4.5-/im Spitzer channel to constrain the CO2 abundance 
could be problematic due to the additi onal contribu- 
tion o f CO opacity at these wavelen gths (jFortnev et al.l 
I2010t IMadhusudhan fc Seagerl I2009T ). For example, al- 
though our lx solar, C/O = 0.88, disequilibrium 
model has a CO2 abundanc e that falls outside the 
IMadhusudhan fc Seagerl (|2009D £ 2 < 1 constraint, the 
model spectrum has a band-averaged 4.5-fim flux that 
falls within 1-sigma of the observed flux, and solar metal- 
licities are still plausible for HD 189733b, at least as 
far as the Spitzer secondary-eclipse data are concerned. 
In fact, Fig. [5] demonstrates that there is not much 
overlap in the solution s that are consistent with the 
IMadhusudhan fc Seagerl d200l C0 2 and H 2 £ 2 < 1 
constraints in terms of metallicity vs. C/O ratio phase 
space, and we consider the H2O abundance to be the 
more reliable constraint, as water has a larger effect on 
the spectrum. 
The bottom line from our modeling is that although 



the Spitzer/IRAC secondary-eclipse observations indi- 
cate that the C/ O ratio is likely less than u nity on 
HD 189733b (see IMadhusudhan fc Seagerl 120091) . atmo- 
spheric models with C/O ratios moderately enhanced 
over the solar value of ^0.55 provide better fits to the 
Spitzer IRAC and MIPS secondary-eclipse data than 
models with solar-like C/O ratios, given our assumed 
nominal thermal profile, unless the atmosphere has a 
subsolar metallicity. More specifically, the flux in the 
longer- wavelength Spitzer bandpasses is very sensitive to 
the water abundance, and the relat i vely l ow H2O abun- 
dance that Madhusudh an fc Seagerl (|2009) infer from the 
secondary-eclipse observations implies low metallicities, 
high C/O ratios, and/or high dayside temperatures (i.e., 
very inefficient heat redistribution) on HP 189733b (see 
also lBarmanll200l IBurrows et aLll2008l) . 

Disequilibrium processes like photochemistry and 
transport-induced quenching do not change this conclu- 
sion because water is efficiently recycled in the middle at- 
mos phere through pho tochemical and kinetics processes 
(see IMoses et al l 120111 for more details), although dis- 
equilibrium chemistry is found to affect the profiles of 
other spectrally important molecules such as CH4, NH3, 
and HCN (see Fig. [3] and iLine et "all 120101 IMoses et all 
1201 ID . We reproduce Spitzer secondary-eclipse photo- 
metric data well for models with metallicities between 
1-5 x solar and C/O ratios between 0.88-1.0. It should 
be kept in mind, however, that our fits are non-unique. 
The C/O ratio that provides the best fit to the obser- 
vations will depend on both the metallicity and thermal 
structure, with higher metallicities and/or lower temper- 
atures requiring higher C/O ratios to remain consistent 
with the data. 

3.1.3. HD 189733b Recent Updates 

The most rec ent HD 189733b Spitz er IRAC secondary 
eclispse data of iKnutson et al.l (|2012f ) imply significantly 
reduced fluxes in the 3.6 and 4.5 /jm channels (see Fig. [6]). 
With this new data, the derived F p i anet / F sta r at 3.6-/zm 
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Fig. 6. — Secondary-eclipse spectra for HD 189733b (in terms of the ratio of the flux of the planet to that of the host star) for models with 
(Left) a C/O ratio of 0.7 and a metallicity of 4 times solar, and (Right) a C/O ratio of 1.01 and a metallicity of 5 times solar. The solid 
black lines represent synthetic spectra from the result s of our disequi librium chemistry models. The red circles with error bars represent the 
3.6 and 4.5 fim Spitzer p hotom etric data from lKnutson eF al. (2012) and the 5.8, 8, 16, and 24 /im data from Knutson ct al. (2007, 2009b) 
and ICh arbonncau ct al. (2008J); the blue circles with error bars represent the older d ata in the 3.6 and 4.5 ^m IRAC channels. The purple 
circles with error bars represent the HST/N ICMOS spectra of Swain ct al. (2009a), and the yellow circles with error bars represent the 
Spitzer IRS spectra of Grillmair ct al. (2008). The green circles without error bars represent our model results convolved over the Spitzer 
instrument response functions. Note that the thermal contrast between the troposphere and stratosphere was reduced in our new models 
to better fit the new IRAC 3.6 and 4.5 fim points, and the resulting synthetic spectra compare better with the IRS spectra than the older 
models shown in Fig. [4] A color version of this figure is available in the online journal. 



no longer exceeds that at 4.5-/im, and in fact the entire 
spectrum more closely resembles a blackbody. If the scat- 
tering hazes that are observed at ultraviolet and visible 
wavelengths (e.g.. iSing etaLl 120111: iHuitson et all 120121) 
are vertically optically thick at infrared wavelengths, 
then the spectrum at secondary eclipse might indeed be 
expected to resemble a blackbody, but a similar result 
could be obtained with optically thin hazes and less ther- 
mal contrast between the stratosphere and troposphere, 
or from a low-metallicity pla net. We have c r eated new 
models to better fit the new iKnutson et al.l (|2012t ) 3.6 
and 4.5 /im photometric points, along with the older 
5.8, 8 16, and 24 urn data from IKnutson et all (|2007l 
2009tJ) and iCharbonneau et al.l (|2008t ). Some new re- 
sults are shown in Fig. H>] When we consider models 
with a reduced AT between the troposphere and strato- 
sphere, the predicted model excess-flux "bulges" at 3.5- 
4 jjtva and 9-12 /im that were apparent in the previous 
models shown in Fig.|3]are correspondingly reduced. The 
new models then comp are better with the IRS spectra of 
IGrillmair et al.l (|2008f) and have the appropriate relative 
fluxes in the 3.6 and 4.5 /Ltm IRAC channels, but the 
predicted amplitude in the near-IR absorption bands is 
smaller than is indi cated by the HST/NICMOS data of 
I Swain et al.l (|2009af ). As with our previous results, we 
find that models with moderately high C/O ratios result 
in the low water abundances that are needed to provide 
better fits to the data than models with solar-like values 
of C/O ~ 0.5; on the other hand, C/O ratios > 1 provide 
too much absorption from methane and HCN at 3.6 and 
8 /itm to be consistent with the IRAC data. Our favored 
models have C/O ratios in the range ~0. 6-1.0. In all 
tests, our models pro vide insufficient absor ption at 4.5 
/urn to account for the IKnutson et al.l ()2012l ) flux in this 
channel, considering the exceeding ly small statist i cal er - 
ror bars cited for the 4.5 fim point. IKnutson et al.l ()2012f ) 
had similiar difficulty with model fits at 4.5 fim. Increas- 
ing the metallicity, with its corresponding increase in the 



CO2 abundance, does not help the situation, as it is still 
difficult to fit the 4.5-/mi point while keeping good fits 
at other wavelengths. On the other hand, the potential 
systematic errors involved with the IRAC analyses, the 
observed variability of the star, and the uncertainties in 
the stellar models needed for calculating F p i anet / F sta r 
for the synthetic spectra all come into play here, and 
fits such as are shown in Fig. [5] may be reasonable when 
all systematic observational and modeling uncertainties 
are considered. In any case, spectral observations by the 
James Webb Space Tel escope or dedicated missions like 
FINESSE and EChO (jSwainl [20ll iTinetti et all [2011 
could help resolve the atmospheric composition and bet- 
ter constrain the C/O ratio on HD 189733b and other 
extrasolar panets. 

3.2. XO-lb Results 

The transiting plane t XO-lb was discovered by 
IMcCullough"et~aLl <l2006h . XO-lb, like HD 189733b, is a 
moderately-irradiated, moderately-hot transiting planet 
that is warm enough that CO should be the dominant 
carbon-bearing constituent in the middle atmosphere of 
the planet. Whether this exoplane t falls into the 01 or 
C2 re gime of iMadhusudhar] (|2012D depends on the C/O 
ratio. Madhusudhan (2012) suggests th at the fit to the 
Sjpitzer/IRAC secondary-eclipse data of Machalek et al.l 
( 2008) is improved for models with C/O > 1, although 
ITinetti et atl pOlOh are able to find good fits to the same 
data (as well as to transit data) for more standard as- 
sumptions about the C/O ratio, depending on the CH4 
line parameters that are adopted. The relative fluxes in 
the 3.6 vs. 4.5-Atm channe ls and 5.8 vs. 8.0-^m channels 
led lMachalek et al.l (|2008D to conclude that XO-lb has a 
thermal inversio n, although the i r mod el fit to the data 
is not very good; [Madhusu dhan! (|2012l ) demonstrates an 
improved fit for m odels with no therm al inversion and 
a C/O ratio > 1. ITinetti et all (|2010D find a good fit 
to the secondary-eclipse data for models with and with- 
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out thermal inversions, for various assumptions about 
species abundances. For the case of no thermal inver- 
sion, the very high flux at 5.8 //m implies a very low 
water abu ndance and/or very high photospheric tem- 
peratures (jTinetti et aJJl2010l : lMadhiisudharJl2012t l. De- 
generacies between the thermal structure and composi- 
tion are cl early plaguing the interpretation of the X O-lb 
data (e.g.. iTinetti et alj|20ld; lMadhusudhanll2012t) , and 
we are unlikely to find a definitive soluti on to the prob- 
lem. We can, however, test whether the iMadhusudhanl 
(2012) high-C/O-ratio scenarios for XO-lb can fit the 
secondary-eclipse data for chemically plausible species 
abundances that take into account disequilibrium pro- 
cesses. 

3.2.1. XO-lb Chemistry 

Figure [7] illustrates how the thermochemical equilibri- 
um composition of XO-lb varies as a function of the at- 
mospheric C/O ratio for our assumed nominal day side 
thermal structure (see Fig.[lJ and an assumed 1.5 x solar 
metallicity. As with HD 189733b, CO is abundant for all 
assumptions of the C/O ratio, whereas H2O and CO2 are 
much more abundant for C/O < 1, and conversely CH4, 
HCN, and C2H2 are muc h more abundant for C/O > 1. 
The spectral analysis of IMadhusudhanl (|2012h suggests 
that the C/O ratio on XO-lb could be greater than that 
onHD 189733b, perhaps even C/O > 1. For a model with 
a thermal structure similar to our adopte d nominal pro- 
file sh own in Fi g. [1] Madhus udhanl (|2012f ) finds good fits 
to the iMachalek et all (|2008H Svitzer IRAC secondary- 
eclipse spectra for the species mole fractions marked by a 
star in Fig. [71 These constraints imply that CO is abun- 
dant, CO2 is much less abundant, and CH4 » H 2 0. A 
C/O ratio near unity can satisfy all the these constraints 
for our adopted thermal profile, except the water abun- 
dance will be over-predicted for a C/O ratio ~ 1. 

Disequilibrium processes do not affect this overall con- 
clusion much, although the vertical profiles of the spec- 
trally active molecules CH 4 , H 2 0, HCN, NH 3 , and C0 2 
are significantly altered by transport-induced quench- 
ing and photochemistry. The results of three of our 
thermo/photochemical kinetics and diffusion models are 
shown in Fig. [7j Further results from one of these dis- 
equilibrium models are shown in Fig. [8] The main re- 
action schemes controlling the abundances in the C/O 
= 0.5 model are the same a s on HD 189 7 33b; a full 
discussion can be found in iMoses et al.l f|2011h and 
iLine et all (|2010f ). The disequilibrium processes affecting 
the molecules do change as C/O exceeds unity, although 
the main reaction schemes controlling the behavior typ- 
ically remain the same. The main differences pertain to 
the relative significance of quenching in controlling the 
mole fractions of different species in the middle and lower 
stratosphere, and the different relative abundances of the 
species at those quench points, which can lead to differ- 
ent "parent" molecules dominating the photochemistry 
in the upper atmosphere. 

For our favored model with 1.5 x solar metallicity, C/O 
= 1.05, and K zz = 10 9 cm 2 s" 1 (Fig. p}, the only 
molecules relatively unaffected by disequilibrium pro- 
cesses (except at very high altitudes) are CO and N 2 . 
Compared with equilibrium, photochemistry strongly de- 
pletes CH4 above ^0.4 mbar and NH3 above ~1 mbar, 
although the NH3 mole fraction in the 1-10 -3 bar region 



is enhanced by t ranspo rt-induced quenching. Scheme (2) 
of IMoses etaLl ()201lD still dominates the CH 4 -» CO 
quenching, but the very high equilibrium CH4 abundance 
at the quench point allows methane to be a more impor- 
tant parent molecule for the subsequent photochemistry, 
with C2H2 production in particular being enhanced in 
the C/O = 1.05 model as compared with the C/O = 0.5 
model. The large H abundance in the upper atmosphere 
drives the conversion of CH4 into C2H2 through schemes 
such as 



2 (H + CH 4 
2 CH 3 + M 
H + C 2 H 6 
C2H5 + M 
H + C2H4 
H + C 2 H 3 



CH 3 + H 2 ) 
C 2 H 6 + M 
C2H5 + H 2 
H + C2H4 - 
C2H3 + H2 
C2H2 + H2 



M 



Net : 2CH 4 + 4H -> C 2 H 2 + 5H 2 , 



(1) 



where M represents any third body. The key to this 
mechanism operating in the CH4 — > C2H2 direction is the 
reaction pair H + CH4 ^> CH 3 + H2 , which can stay bal- 
anced in the lower stratosphere but operates with a net 
imbalance to the right in the upper stratosphere due to 
the large abundance of H atoms from H2O p hotolysis and 
subse quent catalytic destruction of H2 (see IMoses et al.l 

131. 

The greater abundance of CH 4 in the C/O = 1.05 
model also leads to significantly enhanced disequilib- 
rium production of H CN in the ~0 . 1-to-1 0~ 4 bar region 
through scheme (7) of IMoses et al.l f|201 lh . At pressures 
less than 10 -4 bar, the HCN is derived from CO and N2 
through schemes such as 



CO - 
N 2 - 
- 
N + 
C + 

CN - 
- 

OH H 
2(2H 



- hv 

- hv 
H 2 

OH 
NO 
H 2 
H 2 
H 2 
f M 



C + 
2N 
OH - 
NO - 
CN 4 
HCN 
OH - 
H 2 
H 2 + 



H 
H 



+ H 
H 

h H 

M) 



Net : CO + N 2 + 2H 2 -> HCN + H 2 + N + H, 

(2) 

where hv represents an ultraviolet photon. 

Note that photolysis of CO in schemes like the one 
above leads to the photochemical production of H2O 
at high altitudes, despite the low predicted equilibrium 
abundance of H2O in this region of the atmosphere in 
the C/O = 1.05 model. In fact, HCN and H 2 have 
similar mole fractions at high altitudes in the C/O = 
1.05 model because the O from CO photolysis ends up 
largely as H2O, whereas the C ends up in HCN. Both 
H2O and HCN are rec ycled efficiently at high altitudes 
(see IMoses et al.ll201lD . Transport- induced quenching is 
an important disequilibrium process for water in the C/O 
= 1.05 model, acting to smooth out the mole- fraction 
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Fig. 7. — Mixing-ratio profiles for H2O, CH4, CO, CO2, HCN, and C2H2 on XO-lb from our thermochemical-equilibrium models with 
assumed 1.5x solar metallicity but different assumed C/O ratios (ranging from 0.1 to 1.9, incrementing by 0.1 — dashed black lines for 
C/O < 1, solid lines for C/O > 1). The colored lines represent disequilibrium-chemistry results for models with a solar-like C/O ratio of 
0.5 and K zz = 10 9 cm 2 s — 1 (purple), a model with a C/O ratio of 1.05 and K zz = 10 9 cm 2 s — 1 (green), and a model with a C/O ratio of 
1.05 and K zz = 10 8 cm 2 s — 1 (blue). All models have an ass umed atmospheric metallicity of 1.5x solar. The red stars repr esent abundance 
constr aints from a model presented in Madhusudhan (2012) based on the Spitzer/YRAC secondary-eclipse observations of Machalck et al. 
(2008); HCN and C2H2 were not considered in that analysis. A color version of this figure is presented in the online journal. 
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Fig. 8. — Mole-fraction profiles for several important species 
(as labeled) in a thermochemical and photochemical kinetics and 
transport model for XO-lb. Model assumptions include 1.5 X so- 
lar metallicity, a C/0 ratio of 1.05, and K zz = 10 9 cm 2 s — 1 
(see also the green curve in Fig. 0. The dotted lines show the 
thcrmochemical-equilibrium abundances for the species, using the 
same color coding. The overall messiness of this figure indicates 
that disequilibrium chemistry is important on XO-lb. A color ver- 
sion of this figure is available in the online journal. 

profile such that the H2O abundance drops below equi- 
librium values near 0.1-0.01 bar but exceeds equilib- 
rium values above ^10~ 2 mbar. Quenching of water 
occurs through the same overall reaction scheme that 
af fects the quenching of CO and CH4 (e.g., scheme (2) 
in iMoses et al.l 1201 If f . This scheme also operates in the 
C/0 < 1 models, but the quenching of water is more 
obvious in the C/0 > 1 models due to the lower equilib- 
rium H2O abundance. Carbon dioxide is also affected by 
the quenching of CO and H2O, and the increased pho- 
tochemical production of H2O and OH at high altitudes 
leads to a significantly increased CO2 abundance in the 
C/0 = 1.05 model as compared with equilibrium, but 
the overall abundance of CO2 remains much lower than 
that of water. 

As with the equilibrium model, the disequilibrium 
model has m o re H9 O than the spectral analysis of 
iMadhusudhanl (|2012[ ) would indicate, but the resulting 
abundances of CO a nd CH4 are still in line with the 
IMadhusudhanl (|2012D constraints. The disequilibrium 
models underpredict the CO2 abundance on XO-lb, as 
they did for HD 189733b, but again, high CO2 abun- 
dances are not required to fit the thermal-infrared data at 
secondary eclipse (see below). The disequilibrium models 
also indicate that hydrogen cyanide is abundant enough 
that HCN opacity should be included in spectral models. 

3.2.2. XO-lb Spectra 

In Fig. [5J we compare a synthetic spectrum from our 
nomin al disequilibrium mode l with the XO-lb Spitzer 
data of Machalck et al. ( 2008) . The photometric data in- 
dicate a flux at 4.5 /im that is larger than that at 3.6 /im 
and a flux at 5.8 /jm that is larger than that at 8.0 /xm. In 
the absence of a thermal inversion, such flux ratios sug- 
gest that CH4 (an d/or HCN) are more abundant than 
H 2 0, which is why Madhus udhanl (|2012l ) suggested XO- 
lb could be carbon-rich in the first place. Our spectral 
fit for this model with C/O = 1.05 is within 1.7-sigma 
of the data for all the Spitzer bands, with the biggest 
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Fig. 9. — Secondary-eclipse spectra for XO-lb. The solid curves 
represent synthetic spectra from the results of our disequilibrium 
chemistry models with 1.5x solar metallicity, K zz = 10 9 cm 2 s _1 , 
and C/O ratios of either 1.05 (green) or 0.5 (purple). The blue 
circles with error bars represent the Spitzer IRAC observations of 
IMachalek et al.1 JJOOS). The circles without error bars represent 
the model results for C/O = 1.05 (green) or C/O = 0.5 (purple) 
convolved over the Spitzer bandpasses. Note the much better fit 
to the data with the C/O = 1.05 model for this particular thermal 
structure and metallicity. The dashed lines are blackbody curves 
for assumed temperatures of 900 K (bottom), 1250 K (middle), 
and 1600 K (top). A color version of this figure is available in the 
online journal. 

problem being an overestimation of the absorption at 5.8 
/im. Going to a higher C/O ratio, with its correspond- 
ingly lower water abundance, would improve the fit in 
this wavelength band, but the resulting larger CH4 and 
HCN abundances would degrade the fit at 3.6 and 8.0 
/im. Our spectral comparison here demonstrates that 
carbon-rich models with plausible species abundances 
from disequilibrium-chemistry predictions can be consis- 
tent with the Sp itzer secondary - eclips e data, validating 
the suggestion of IMadhusudhanl (|2012l ). 

Given that XO-lb is cool enough that TiO is not 
expected to be present in significant quantities in the 
stratosphere, even for C/O ratios < 1 (and thus TiO 
would not be available to produce a stratospheric ther- 
mal in version; iHubenv et ail 120031 iFortnev et alll2006bl 
2008a), the high C/O-ratio, non-inverted scenario for 
XO-lb has some theoretical advantages. However, we 
must again emphasize that our models are not unique. 
The results will be sensitive to the assumed metallic- 
ity and thermal structure, for example. We note that 
the relative 3.6-/zm/4.5-/zm and 5.8-/im/8.0-/im flux ra- 
tios from the XO-lb secondary-eclipse data cannot be 
explained with solar-like C/O ratios unless temperatures 
at ~0.1-l-bar are much hotter (>2000 K) than is sug- 
gested by theoretical predictions or u nless the planet 
has a stratospheric therma l inversion (cf. IMachalek et al.1 
l2008t ITinetti et aL|[20ToT) . In fact, the observed 5.8- 
/im/8.0-/im flux ratio is difficult to explain with any 
model — whether thermally inverted, anomalously hot 
at depth, or possess i ng a high C/ O ra t io (see Fig. El 
Mach alek et al.ll2008l ITinetti et aLll2010L IMadhusudhanl 
I2012h — and additional high-precision photometric or 
spectral data from ground-based or space-based obser- 
vations would be highly desirable in helping constrain 
chemical and physical properties of this unusual and in- 
triguing exoplanet. 
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3.3. WASP-121 Results 

The tran s iting planet WASP-12b, discovered by 
iHebb et al] (|2009l ). is much hotter than either HD 
189733b or XO-lb (see Fig.[TJ). As one of the most highly 
irradiated hot Jupiters discovered to date, WASP-12b 
would fall into the upper part of e ither the 02 or C2 
regimes of the iMadhusudhanl (j2012T ) hot-Jupiter classi- 
fication scheme, depending on the C/O ratio. WASP- 
12b is hot enough that titanium would not be cold- 
trapped into conde nsed calcium-titanium oxides (e.g., 
iLodderd 12002. 2010]) or other condensed phases. Thus, 
for C/O ratios less than unity, TiO vapor would be a 
major titanium-bearing phase that would be expected 
to survive into the stratosphere of WASP-12b, where 
it could absorb at optical wavelengths, heat the at- 
mosp here, and create a s t ratospheric thermal inversion 
(e.g.. iHubenv et al l [200l iFortnev et all I2006U l2008at 
iBurrows et al . 2008; S piegel et al . 2009) such as has been 
suggested for HD 209458b, HAT-P-7b, and several other 
transiting planets based on Syitzer observations (e.g. , 
Burrows et all 120071 120081: iKnutson et al.l [200l | 2009a|: 
Swain et al.ll2009bt|Madhusudhan fc Seagerl|2009l l2010at 
Spiegel fc Burrowsll2010l; iChristiansen et al.ll2010h . 
However, Madh usudhan et al.l (I2011afl conclude that 



a strong stratospheric thermal inversion in the 0.01-1 
bar region of WASP- 12b can be statistically ruled out 
due to a poor fit to the seco ndary-eclipse Spitzer/IRAC 
data of iCampo et al.l (|2011) and ground-based data o f 
iCroll et al.l (|201lD . Instead. iMadhusudhan et all (|2011aft 
suggest that the atmosphere of WASP-12b is enriched 
in carbon relative to solar elemental abundances, with 
WASP- 12b spectra exhibiting strong absorption due to 
CO and CH4, and much weaker absorption due to H2O, 
as would be expected from an atmosphere with a C/O 
ratio greate r than unity. In fact , the mo del-data com- 
parisons of IMadhusudhan et all (|2011aD suggest that 
a solar-like ratio of C/O = 0.54 for WASP-12b can 
be ruled out statistically at the 4.2cr-significance level, 
whereas chemical-equilibrium models assuming C/O = 1 
yield abundances of CO, H 2 0, and CH4 that are con- 
sistent with s e conda ry-eclipse observations. Recently, 
IMadhusudhanl (|2012f) updated the WASP-12b secondary- 
eclipse analysis to include HCN and C2H2 in the spectral 
modeling, and their conclusions with respect to the high 
C/O ratio have not changed. Moreover, HS T Wide Field 
Came ra 3 (WF3) spectral observations by I Swain et al.l 
(2012) confirm the suggestion of a very low water abun- 
dance on WASP-12b, as well as a low TiO abundance, 
furtherin g the infere nce of a hi gh C/O ratio f or this 
planet. ISwain et all ((2012) and IMadhusudhanl (|2012l ) 
demonstrate that the WFC3 secondary-eclipse observa- 
tions are well reproduced by high C/O ratio m odels. 
However, it is worth noting that ISwain et al.l (j2012|) favor 
a very-low- metallicity solution over a high-C/O-ratio so- 
lution in explaining the WFC3 observations because the 
high C/O ratio models provide a worse fit to the tran- 
sit observations, especially at the short-wavelength end 
near 1.1-1.3 /zm. The degree to which Rayleigh scatter- 
ing could be affecting the shorter wavelengths was not 
discussed. 

In any case, the IMadhusudhan et ail (|2011aD analy- 
sis provided the first strong evide nce for a high C/O 
ratio on a transiting exoplanet. [Madhusudhan et al.l 



(|201 lbf) further emphasized that TiO would not be abun- 
dant in giant-exoplanet atmospheres where C/O > 1, 
thus providing a theoretically consistent explanation for 
the la ck of a stratosphe r ic the rmal inversion on WASP- 
12b. IKopparapu et al.l ()2012D followed up this analy- 
sis by developing photochemical models for WASP- 12b 
for two different assumptions of the atmospheric C/O 
ratio; the molecules HCN and C2H2 were included in 
these models based on a suggestion bv lMoses fc Visscherl 
(|2011D that these sp ecies would b e important for high 
C/O r at ios (see alsolLo ddcrs 2010; Madhusudha n et al.l 
I2011bf) . IKopparapu et all (j2012( ) found that C 2 H 2 and 
HCN were even more abundant than CH4 in their pho- 
tochemical model with C/O = 1.08, and they sug- 
gested that C2H2 rather than CH4 is responsible for the 
strong 8-/im abso r ption on WASP-12b. We expand the 
IKopparapu et a l. (2012) modeling by investigating the 
equilibrium and disequilibrium composition of WASP- 
12b for a variety of C/O ratios and by calculating the 
spectral consequences of abundant HCN and C2H2 for 
the carbon-rich cases. 

3.3.1. WASP-12b Chemistry 

Figure [TU] illustrates how the thermochemical equili- 
brium composition of WASP-12b varies as a function of 
the atmospheric C/O ratio for our assumed nominal day- 
side thermal structure (see Fig. Q} and an assumed 3x 
solar metallicity. Because of WASP-12b's high tempera- 
tures, the chemical-equilibrium abundance of methane in 
the photosphere of WASP- 12b is predicted to be smaller 
than that on the cooler HD 189733b or XO-lb, for similar 
assumptions of atmospheric metallicity and C/O ratio. 
Methane, with its relatively weak C-H bond, is less stable 
at high temperatures and low pressures than CO, which 
has a very strong carbon-oxygen bond. As with the other 
planets we have investigated, CO is again abundant for 
all assumptions of the C/O ratio, whereas H2O and CO2 
are only important for C/O < 1, and HCN and C2H2 be- 
come more important for C/O > 1. At the high tempera- 
tures predicted for WASP-12b, HCN — which also has a 
strong carbon-nitrogen bond — quickly replaces CH4 as 
the second-most abundant carbon-bearing molecule (be- 
hind CO) as the pressure decreases, even for C/O < 1. 
With increasing C/O ratio, C2H2 also becomes a major 
equilibrium constituent whose abundance can eventually 
exceed that of CO when C/O ^ 1. The implied large 
methane a bundance required to fit the s econdary-eclipse 
data in the IMa dhusudh an et all (|2011af) statistical anal- 
ysis suggests a C/O ratio of ~2 or larger for WASP- 
12b (see Fig. ITU1) ; however, this conclusion changes when 
HCN and C2H2 opacities are included in the spectral 
models. For example, as is discussed in Section 13.3.21 at 
least part of the opacity hitherto attributed to CH4 on 
WASP-12b could instead result from HCN and/or C 2 H 2 , 
which are abundant species in carbon-rich atmospheres 
with C/O ratios > 1 (see also IKopparapu et all 120121 ; 
IMadhusudhanl l20Tl) . 

Disequilibrium processes like photochemistry and 
transport-induced quenching result in relatively minor 
changes in the profiles of the dominant species in our 
C/O = 1.015 model, except for HCN, C2H2, and atomic 
species (see Fig. [11]). Hydrogen cyanide becomes a par- 
ent molecule of other photochemical products in much 
the same way as NH3 and CH4 do in cooler atmospheres. 
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Fig. 10. — Mixing-ratio profiles for H2O, CH4, CO, CO2, HCN, and C2H2 on WASP-12b from our thermochemical-equilibrium models 
with assumed 3x solar metallicity but different assumed C/O ratios (ranging from 0.1 to 1.9, incrementing by 0.1 — dashed black lines 
for C/O < 1, solid lines for C/O > 1). The colored lines represent disequilibrium-chemistry results for models with 3x solar metallicity 
and a solar-like C/O ratio of 0.5 (purple), a model with 3x solar metallicity and a C/O ratio of 1.015 (green), and a model with 5x solar 
metallicity and a C/O ratio of 2.0 (blue). All models have an assumed vigorous K zz = 10 10 cm 2 s _1 , which might b e expected for such a 
highly irradiated planet. The red horizontal line segments represent abundance constraints from models presented in Madhusudhan ct al. 
I|2011afl based on the Spitzer /IR AC secondary-ecli pse observations of [Campo ct al. (201l|) and the ground-base d Ks, H, and X-band 
secondary-eclipse observations of Croll et al. (2011). Note that HCN and C2H2 were not considered in the [Madhusudhan c t al.l l|2011aA 
analysis, and the red star in the HCN figure is from the present analysis and indicates the HCN abundance that best reproduces the 
secondary-eclipse data for our adopted nominal thermal structure shown in Fig. [T] A color version of this figure is presented in the online 
journal. 
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Fig. 11. — Mole-fraction profiles for several important species 
(as labeled) in a thermochemical and photochemical kinetics and 
transport model for WASP-12b. Model assumptions include 3x 
solar metallicity, a C/O ratio of 1.01 5, and K zz = 10 10 cm 2 s _1 
(see also the green curves in Fig. 1101 . The dotted lines show the 
thermochemical-equilibrium abundances for the species, using the 
same color coding. Note that both HCN and C2H2 are more abun- 
dant than CH4 throughout the middle atmosphere. A color version 
of this figure is available in the online journal. 

Due to the high temperatures and large overall H abun- 
dance at high altitudes, kinetic loss of HCN via H + HCN 
— > CN + H2 dominates over HCN photolysis as the main 
photochemical loss process. The CN radicals produced 
can recycle the HCN or become photolyzed themselves, 
leading to C + N. Atomic N can then react either with 
NO to form N 2 + O or with CN to form N 2 + C, and the 
atomic carbon can react further to eventually produce 
C2H2. The dominant C2H2 photochemical production 
scheme in our C/O = 1.015 model is 



H 2 + M -> 2 H 
2 (H + HCN -> CN 
2 (CN + hp -> C + 
2 (C + H 2 CH 



CH + H 2 
3 CH2 + H2 
CH + CH 3 
H + C2H3 
N + H 2 
NH 
N -I 
3(2H 



+ 
NO 



d CH 2 
CH 3 
C2H3 
C2H2 
NH -} 
NO -f 
N 2 + 



- M 

- H 2 ) 
N) 
-H) 

+ H 
r- H 
+ H 
+ H 2 
H 
H 




M — > H 2 + M ) 



Net : 2 HCN C 2 H 2 + N 2 . 

(3) 

The acetylene abundance increases with altitude in the 
^2-mbar to 2-/ibar region due to this HCN photochem- 
istry, while the HCN abundance decreases. The nitrogen 
from HCN loss largely ends up as N 2 at pressures greater 
than ~1 /ibar and as atomic N at pressures less than ~1 
/xbar. In fact, no molecules (including H 2 ) survive at alti- 
tudes much above ~1 /^bar, as the intense ultraviolet flux 
and large H abundance attack most molecular bonds and 
favor atomic species. Some photochemical production of 
H 2 0, CH 4 , and NH 3 occurs from the CO and HCN pho- 
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Fig. 12. — Secondary-eclipse spectra for WASP-12b. The solid 
lines represent synthetic spectra from the results of our disequi- 
librium chemistry models with 3x solar metallicity, K zz = 10 10 
cm 2 s , and C/O ratios of either 1.015 (green) or 0.5 (purple). 
The blue circles w i th err or bars represent the Spitzer / IRAC data 
from Campo ct al. (2011) and ground-based data from Croll ct al. 
H201H 1. The circles without error bars represent the model results 
for C/O = 1.015 (green) or 0.5 (purple) convolved over the Spitzer 
bandpasses. Note the better fit with the C/O = 1.015 model for 
this assumed metallicity and thermal structure, particularly at 2.1, 
3.6, and 5.8 (im. The dashed curves represent blackbodies with 
temperatures of 1800 K (bottom), 2500 K (middle), and 3000 K 
(top). A color version of this figure is available in the online jour- 
nal. 

tochemistry, but these species remain below ppm levels 
at pressures less than ~0.1 bar and only become signif- 
icant at high pressures and in a narrow region centered 
near ~1 /ibar where photolysis of CO and N2 is signifi- 
cant. In all, disequilibrium processes have less of an effect 
on hotter planets like WASP- 12b than they do on cooler 
planets like HP 189733b and XO-lb (see a lso Lian g et all 
2001 [200l IZahnle et alj|2009bt ILine et al.| [2010l |2011at 
Moses et alj|2011t iVenot et al.ll2012D . 

Our photochemical models for a high C/O ra- 
tio on WASP- 1 2 b dif fer somewhat from those of 
IKopparapu et all (|2012l ). Their upper atmosphere is 
cooler, which tends to reduce reaction rates. Transport- 
induced quenching of CH 4 and H 2 therefore become 
correspondingly more important, as expected, leading to 
increased abundances of H2O and CH4. Their adopted 
temperature-pressure profile and C/O ratio of 1.08 also 
result in C2H2 being slightly more abundant than HCN 
in equilibrium throughout the photosphere, which can 
affect the photoche mistry. However, the HC N and C2H2 
vertical profiles of IKopparapu et all (|2012t ) are almost 
completely unaffected by photochemistry, in contrast to 
our models, for reasons that are unclear but likely have to 
do with different reaction mechanisms. Their main con- 
clusion that HCN and C2H2 are important on WASP-12b 
for the case of C/O > 1 under both equilibrium and dis- 
equilibrium condit ions is qualitatively supported by our 
models (see alsolLodderdl20l0HM oscs fc Visscherl l20ll 
iMadhusudhan et aLll201lU lMadhusudhanll2012fh 



3.3.2. WASP-12b Spectra 

Like methane, both HCN and C2H2 have bands that 
fall within the Spitzer/IRAC 8.0-/mi spectral bandpass, 
and to a lesser extent within the short- wavelength end 
of th e 3.6-/xm spectral bandpass (e.g., iShabram et al.l 
l2011f ). We find that all three molecules contribute to 



18 



Moses et al. 



the calculated flux in the 3.6 and 8-/xm channels, but for 
our nominal thermal structure and C/O ratio ~ 1, the 
channel-integrated flux at these wavelengths is particu- 
larly sensitive to the HCN mole fraction. For a 3x solar 
metallicity and our adopted thermal structure, C/O ra- 
tios just over unity result in an HCN abundance that 
provides a good fit to the Spitzer/IRAC 3.6-/im band 
depth, while providing slightly too much absorption in 
the lower-signal-to-noise 8.0-/xm band. This result will 
be model-dependent, particularly in terms of the relative 
contributions of different molecules for different thermal 
structures, but Fig. ll2l demonstrates that our model with 
C/O = 1.015 provides a significantly better fit (with £ 2 
< 2.6) to the Spitzer data than the model with C/O 
= 0.5 (with £ 2 < 10.2); this resu lt is consistent with 
the det ailed statistical analysis of lMadhusudha n~et al.l 
(|2011af) . The C/O = 1.015 model has much less ab- 
sorption at 5.8 /xm and much more absorption at 3.6 
/jm than the C/O = 0.5 model, and the predicted rela- 
tive strengths in all four Spitzer/IRAC bands are much 
more consistent with observations for the higher C/O 
ra tio model. We ther e fore co ncur with the suggestion 
of iMadhusudhan et al.l (|2011aD that WASP-12b could be 
a carbon-enriched planet, although in our model it is 
HCN rather than CH4 that provides the dominant opac- 
ity source in the 3.6 and 8.0 /xm channels (see also Sec- 
tion 4 and lMadhusudhanll2qi2l ). Note that a high C/O 
ratio model for WASP-12b is also more consistent with 
th e recent HST WF C3 secondary-eclipse observations 
of iSwain et al.l (]2012l ) in which a very low water abun- 
dance is needed to explain the nearly featureless 1-0-1-7 
micron spectrum of WASP-12b (e.g., ISwain et al.l 120121 : 
iMadhusudhanllMl . 

The large predicted HCN abundance on a putative 
carbon- rich WASP- 12b is also l ikely to affec t its n ear- 
infrared transit spectrum, where ICowan et all (|2012[ ) ob- 
serve a larger transit depth in t he 3.6-/xm Sp i tzer c han- 
nel than in the 4.5 /xm channel. ICowan et al.l ()2012D did 
not consider the possible influence of HCN on the 3.6- 
/im absorption, and in fact only considered enhanced CO 
abundances as a consequence of a potential carbon-rich 
scenario. Future transit simulations should consider the 
influence of HCN and C2H2. At 1-2 /xm, the spectrum 
in the high C/O ratio case is much more featureless than 
the solar C/O ratio case, where water absorption bands 
are prominent. 

3.4. CoRoT-2b Results 

The trans i ting planet CoRoT-2b, discovered by 
lAlonso et al.l (|2008lh has often been call ed a "misfit" 
because of its unusually large radius (e.g.. lAlonso et al.l 
[2001 IGillon et al.1[20lfA IGuillot fc Hlrvelll2011l) and be- 
cause traditional models that assume solar elemen- 
tal abundances cannot explain its unusual flux ratios 
from secondary-eclipse observations in the Spitzer/IRAC 
channels at 3.6, 4.5, and 8.0 /xm, regardless of whether 
the planet is ass umed to have a stratospheric thermal in - 
version or not (|Gillon et al.ll2010b IDeming et al.l 120111 ). 
The main problem is an anomalously low flux at 8.0 
zxm and an unusuall y large 4.5-/xm/8.0-/xm flux ratio. 
IMadhusudhan! (|2012l ) suggests that all the secondary- 
eclipse observations for CoRoT-2b can be reproduced 
to within their 1-a unce rtainties for chemicall y plausible 
models with C/O > 1; IMadhusudhan! (|2012|) also finds 



good solutions for oxygen-rich models, but those mod- 
els tend to have implausibly large CH4/CO abundance 
ratios. Because there is no Spitzer/IRAC data available 
for CoRoT-2b from the 5.8-/xm IRAC channel, model so- 
lutions foj_CoRo^^b will tend to be more degenerate, 
but if IMadhusudhan! (|2012t ) is correct in his suggestion 
of a carbon-rich atmosphere, then CoRoT-2b would fall 
within the C2 regime in their classification scheme. 

3.4.1. CoRoT-2b Chemistry 

Figure [13] illustrates how the thermochemical- 
equilibrium composition of CoRoT-2b varies as a func- 
tion of the atmospheric C/O ratio for our assumed nom- 
inal dayside thermal structure (see Fig. [T]) and an as- 
sumed 0.5 x solar metallicity. The low metallicity we 
are considering here is required to prevent too much ab- 
sorption from CO in the Spitzer/IRAC channel at 4.5 
/xm for our assumed thermal structure; a warmer atmo- 
sphere would allow higher metallicities while still being 
consistent with the data. As with all the other planets we 
have considered, CO is abundant in chemical equilibrium 
for all assumptions of the C/O ratio, whereas H2O and 
C0 2 are only abundant for C/O < 1, and CH 4 , HCN, 
and C2H2 become more abundant for C/O > 1. For a 
model with a thermal structure s imilar to our adopted 
nominal profile show n in Fig. [H IMadhusudhan! (I2012T) 
finds good fits to the lAlonso et al.l (120101V IGillon et al.1 
( 2010), and IDeming et al.l ()2011[ ) secondary-eclipse data 
for the species mole fractions marked by a star in Fig. 1131 
These constraints imply that CO and CH4 are abundant, 
whereas H2O and CO2 remain below ppm levels. It ap- 
pears that a C/O ratio near 1.1 can readily satisfy all 
the these constraints. The relatively low upper atmo- 
spheric temperatures adopted for our nominal model sug- 
gest that photochemistry will be important on C0R0T- 
2b, as we found for XO-lb and HD 189733b. 

Indeed, Fig. [14] shows that disequilibrium chemistry 
is important for our CoRoT-2b model with an assumed 
C/O ratio of 1.1, a subsolar metallicity of 0.5x solar, 
and K zz = 10 9 cm 2 s . Transport-induced quenching 
of water is very important, enhancing the H2O photo- 
spheric column abundance by about an order of magni- 
tude, and CO photochemistry further enhances the H 2 
abundance in the upper atmosphere by many orders of 
magnitude. Photochemical production of water occurs 
via CO photolysis in the upper stratosphere to produce 
C + O, followed by reactions of O with H2 to form OH, 
which reacts with H2 to form water. The atomic C can 
react with NO to form CN and eventually HCN (see 
scheme (2) above) or the C can react with H2 to form 
CH, which eventually goes on to produce C2H2. We are 
even seeing a non-trivial formation rate of small C^II X 
compounds at high altitudes in the model from reaction 
of C and/or CH with C2H2; these G^H X species could be 
soot precursors, although benzene and complex nitriles 
themselves never exceed ppb levels in this constantly- 
illuminated model. Transport-induced quenching also af- 
fects CH4, HCN, NH3, CO2, and many other species in 
the middle atmosphere, whereas photochemistry affects 
profiles in the upper atmosphere. Although photochem- 
ical production of CO2 greatly enhances its abundance 
in the upper atmosphere, the column abundance in the 
photosphere is still quite small, and CO2 should have lit- 
tle influence on spectra. The chemistry of a carbon-rich 
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Fig. 13. — Mixing-ratio profiles for H2O, CH4, CO, CO2, HCN, and C2H2 on CoRoT-2b from our thermochcmical-equilibrium models 
with assumed 0.5x solar metallicity but different assumed C/O ratios (ranging from 0.1 to 1.9, incrementing by 0.1 — dashed black lines 
for C/O < 1, solid lines for C/O > 1). The colored lines represent disequilibrium-chemistry results for models with a solar-like C/O ratio 
of 0.5 and K zz = 10 9 cm 2 s _1 (purple), a model with a C/O ratio of 1.1 and K zz = 10 9 cm 2 s _1 (green), and a model with a C/O ratio of 
1.1 and K zz = 10 s cm 2 s _1 (blue). All models have an as sumed atmospheric metallicity of 0.5x solar. The red stars r epresent abundance 
con straints from a prelim inary model ( Madhusudhan 2012) based on the Spitzer/IRAC seco ndary-eclipse observa tions ofjGillon ct al.| 1 120101 ) 
and Doming: ct al. (2011) and the ground-based secondary-eclipse Ks-band observations of Alonso ct al. (2010), in which HCN and C2H2 
were not considered in the analysis. A color version of this figure is presented in the online journal. 
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Mole Fraction 

Fig. 14. — Mole-fraction profiles for several important species 
(as labeled) in a thermochemical and photochemical kinetics and 
transport model for CoRoT-2b. Model assumptions include 0.5 X 
solar metallicity, a C/O ratio of 1.1 , and K zz = 10 9 cm 2 s _1 
(see also the green curves in Fig. I13D ■ The dotted lines show the 
thermochemical-equilibrium abundances for some of the species, 
using the same color coding. Note the enhanced disequilibrium 
abundances of HCN, H2O, NH3, and CO2 over equilibrium values. 
A color version of this figure is available in the online journal. 
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Fig. 15. — Secondary-eclipse spectra for CoRoT-2b. The solid 
lines represent synthetic spectra from the results of our disequi- 
librium chemistry models with 0.5x solar metallicity, K zz = 10 9 
cm 2 s _1 , and C/O ratios of either 1.1 (green) or 0.5 (purple). 
The blue circles w i th err or ba rs represent th e Spitz er/IRAC data 
from Gillon ct al. (2010) and Doming et al. (2011) and ground- 
based date~^om~^^ns^er|al]l[20^ The circles without error 
bars represent the model results for C/O = 1.1 (green) or 0.5 (pur- 
ple) convolved over the spectral bandpasses. Note the better fit 
with the C/O = 1.1 model at 3.6 fim for this particular assumed 
metallicity and thermal structure. A color version of this figure is 
available in the online journal. 

CoRoT-2b is similar to that of HD 189733b a nd XO - 
lb, which is covered in more detail in Sections 13.1.11 & 
13.2.1) the chemistry of a solar-like-composition C0R0T- 
2b i s similar to that of othe r "warm" transit ing planets 
(see iMoses et all 120111 and iLine et al.l 120101 for further 
details) . 

3.4.2. CoRoT-2b Spectra 

Figure [15] shows how two models with diffe rent C/O ra- 
tios co mpar e with the Spit z er/IR AC data of lGillon et al.l 
(|2010D and iDeming et al.l ([201 If) and the ground-based 
data of lAlonso et al.l (|2010D at secondary eclipse. Both 



models exhibit too much absorption at 4.5 /im, which 
demonstrates why we did not investigate models with 
higher metallicities for this particular thermal struc- 
ture. Both models exhibit less flux in the 8-fim chan- 
nel than in the 4.5 /zm channel, in contrast to obser- 
vations, which further emphasizes the un usual observed 
4.5-/xm/8-/xm flux ratio for this planet (see IDeming et al.l 
1201 It IMadhusudhanl 120121) . although both models are 
within ~1.5-cr of the uncertainties at these wavelengths. 
The main difference between the models lies in the 3.6- 
/im/4.5-/im flux ratio and in the overall fit in the 3.6- 
fim channel, for which the C/O = 1.1 model provides a 
significantly better fit to the existing secondary-eclipse 
data. In all, this carbon-rich model provides a statis- 
tically better fit (£ 2 = 1.3) than the solar-composition 
model s pre sented here (£ 2 = 7.2 ) and in iMadhusudhanl 
(|2012D and IDeming et all (|201lD . while still having the- 
oretica lly consistent species abundances. The sugges- 
tion of IMadhusudhanl (|2012D that CoRoT-2b could have 
a carbon-rich atmosphere is therefore viable, although 
further observations at additional wavelengths would be 
required to fully test the suggestion. 

As an aside, we note that unlike the case for the very 
hot WASP-12b or the much cooler XO-lb, our nomi- 
nal thermal profile for CoRoT-2b first crosses the en- 
statite and forsterite condensation curves right within 
the planet's photosphere (for solar-like elemental com- 
positions) , which suggests that silicate clouds could have 
more spectral consequences on CoRoT-2b than they do 
on much hotter or cooler planets. Opacity from thick 
clouds can reduce the depth (and broadness) of absorp- 
tion bands and make the observed planetary spectrum 
appear more like a blackbody. We note, however, that 
the observed secondary-eclipse data for CoRoT-2b do 
not much resemble a blackbody, particularly in terms 
of the deep absorption in the 8-/im channel. If clouds 
were present near ~0.1 bar or deeper, they would have 
relatively little effect on the photochemistry of C, N, 
and O-bearing species. High-altitude clouds located at 
10~ 3 mbar or above, on the other hand, could shield 
species from photolysis, with some interesting conse- 
quences. We also note that the host star for CoRoT-2b 
is very young and spe ctrally active (jBouchv et al.l [20081 : 
Guillot & Havel 2011), which would result in a high x- 
ray and EUV flux and enhanced charged particle fluxes 
that would be conducive for ion chemistr y in the at- 
mosph ere of CoRoT-2b. As mentioned by IMoses et al.l 
(|2011l) . Titan-like ion chemistry at high altitudes could 
enhance the destruction of N2 (and perhaps CO) and 
enhance the production of complex hydrocarbons, which 
could further influence spectral behavior. 

4. IMPLICATIONS WITH RESPECT TO SPECTRA 

The atmospheric composition of carbon-rich planets 
differs considerably from that of solar-composition plan- 
ets, and those differences have spectral implications, as 
is obvious from Figs. [4] [6] [H [T2l and [15J Photochem- 
ical production of HCN and C2H2 is efficient enough 
on cooler exoplanets (even those that are only mildly 
carbon-enhanced) that these species could affect spectra 
in certain wavelength regions. Acetylene would provide 
opacity mainly in the 13-15 /jm region, and to a lesser ex- 
tent in the ~3 and 7-8 /im regions (with the latter bands 
falling within the Spitzer 3.6 and 8.0-/im channels). Hy- 
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drogen cyanide would provide opacity mainly in the ^3- 
/im, 6.5-7.8 /im, and 12-16 fxm regions (with broader 
wavelength ranges for higher temperatures). Even on 
the hottest transiting exoplanets, species like HCN can 
have a major influence on the total atmospheric opacity 
if the atmospheric C/O ratio exceeds unity. For exam- 
ple, Fig. [TBI shows the major influence that HCN has on 
the infrared spectra of our nominal carbon-rich WASP- 
12b model with C/O = 1.015. Although C 2 H 2 has rel- 
atively little influence on the spectrum, only becoming 
noticeable in the ~3-/im region, HCN opacity completely 
dominates in several wavelength regions, even surpassing 
the contribution of CH4, CO, H2O, and other molecules 
at ^2.4-4.1 /im and in the 6-10 /im region. Note that 
CH4 opacities at high temperatures are poorly known 
and may be underestimated by as much as two orders 
of magnitude or more in our spectral calculations. If we 
approximate this effect by increasing the CH4 abundance 
by a factor of 100, we see that CH4 opacity could still 
be important on WASP-12b, but that does not change 
the fact that HCN is likely a major opacity source in 
the Spitzer 3.6 and 8.0 /im channels. Clearly, HCN is a 
potentially important molecule that can contribute to ex- 
oplanet spectral behavior and should be included in anal- 
yses of exoplanet spectra; C2H2 and NH3 are important 
disequilibrium species under some conditions and should 
also be considered if those conditions are met (e.g., po- 
tentially high C/O ratio or cooler atmosphere; see also 
iMoses et alll20lU iKopparapu et al.ll2012| ). 

Note, however, that to truly determine relative 
molecular abundances from exoplanet spectra, more 
detailed information on the line parameters (includ- 
ing intensities) from hot bands is needed in the 
molecula r spectroscopy databases like GEISA and HI- 
TRA N (|Jacquinet-Husson et all 120081 : iRothman et al.l 
2009). There is currently insufficient published high- 
temperature information on CH4, HCN, and C2H2 to re- 
liably evaluate their potential quantitative contributions 
to the spectra. 

Carbon dioxide can influence the spectra, particu- 
larly near ~2 /im, ~4.3 /im, and in the 14-16 /jm re- 
gion. Our carbon-enhanced nominal disequilibrium mod- 
els that provide good fits to the Spitzer/IRAC data 
tend to have CO2 abu ndances lower than those derive d 
from the analyses of |Madhusudhan & Scagcr (2009), 
iMadhusudhan et al.l (|2011al ). and Madhusudhanl (|2012f l 
for these four exoplanets (see Figs. EJ [H pi & [TJ) . Of all 
the Spitzer/IRAC channels, carbon dioxide has the most 
influence on the flux in the 4.5-/im channel, and Figs. 01 
[H O & [JJ] demonstrate that our band-integrated fluxes 
at 4.5 /im do indeed fall slightly above the IRAC-derived 
fluxes at this wavelength for HD 189733b, XO-lb, and 
WASP-12b (CoRoT-2b being the exception). Because 
other molecules such as CO contribute to the flux in this 
IRAC channel, an improved fit would be obtained for 
slight increases in the assumed metallicity of these plan- 
ets; therefore, our low derived CO2 abundances do not 
seem to be a major problem for the Spitzer IRAC fits. 

That is not the case at near-IR wavelengths, however, 
at least for HD 189733b. Note from Figs. H [H &H 
that our predicted CO2 mole fraction for all four of the 
exoplanets investigated never exceeds ppm levels, regard- 
less of the C/O ratio, for both equilibrium and disequilib- 
rium chemistry. Even for 10 x solar metallicities, Fig. 1171 
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Fig. 16. — The contribution of various molecules in our 3x so- 
lar metallicity, C/O = 1.015 disequilibrium-chemistry model to 
the predicted secondary-eclipse spectra for WASP-12b. The green 
curve rep re sen ts synthetic spectra from our nominal model (see 
Figs. 1111 fc E) , whereas the blue curve is the same model with 
C2H2 eliminated. The two curves overlap almost exactly, illus- 
trating the relative insignificance of C2H2 as an opacity source in 
this model. For the magenta curve, the contribution of HCN has 
been eliminated from the nominal model; a comparison of the ma- 
genta curve with the blue and/or green curves illustrates the dom- 
inance of HCN opacity in controlling the spectra behavior in many 
wavelength regions in our model. The black curve is a synthetic 
spectrum without the contribution of HCN, but in which the CH4 
abundance has been increased by a factor of 100. The red circles 
with error bars represent the Spitze r and ground- based seco ndary- 
eclipse data of ICampo et al.l <2011f l and ICroII et al.l 1120111 1. The 
dashed lines are blackbody curves for assumed temperatures of 
1800 K (bottom), 2500 K (middle), and 3000 K (top). A color 
version of this figure is available in the online journal. 

demonstrates that carbon dioxide is not expected to be a 
major product on hot Jupiters, and photochemistry does 
not change this conclusion. 

The relative unimportance of CO2 in our chemical 
models is in stark contrast to the conclusions based 
on secondary-eclipse observations of HD 18 9733b in the 
near- infrared obtained with HST/NICMOS (|Swain et al.l 
l2009al) . The constraints for CO2 based on this 
datase t range from a de rived mole fraction of 10~ 6 - 
10" 1 (jSwain et al.l [2009al). to a best fit of ~7xl0~ 4 
Madh usudhan fc Seagerll2009D . to a best fit of -2 x 10~ 3 
Lee et al.l|2012D. to derived values of (1.7-6.7) x 10" 3 
Line et al.1 1201 lbT ). Such large CO2 abundances are 
inconsistent with the Spitzer/IRAC data in the 4.5- 
/zm channel, whic h has led to the suggestion o f atmo- 
spheric variability (Ma dhusudhan fc Seager 2009) . How- 
ever, note that the retrie vals from this da ta set suggest 
CQ 2 /H 2 ratios of -3 (jLee et al.l I2012D or even ~30 
(jLine et al.l 1201 lO l . which is not likely for a hydrogen- 
dominated planet. In fact, to get CO2 abundances 
comparable to those of H2O under thermochemical- 
equilibrium conditions for our nominal thermal struc- 
ture, HD 189733b would have to have a metallicity of 
order ^2000, at which point H2 would no longer be the 
dominant constituent by number — a suggestion seem- 
ingly inconsistent with t he planet's overall density and 
mass-radius relationship (jMarlev et al.ll2007T ). 

Disequilibrium chemistry cannot help resolve this 
problem. The dominant production mechanism for CO2 
is CO + OH — > CO2 + H, which is indeed effective on 
the dayside of HD 189733b due to OH production from 
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Fig. 17. — The variation in the thermochcmical-cquilibrium abundances of key spectrally active species as a function of temperature at 
0.1 bar, which is a typical pressure predicted to be within the infrared photospheres of hot Jupiters, for lx solar metallicity and C/O = 
0.5 (top left), for lOx solar metallicity and C/O = 0.5 (top right), for lx solar metallicity and C/O = 1.1 (bottom left), and for lOx solar 
metallicity and C/O = 1.1 (bottom right). Note the relatively low CO2 abundance in relation to H2O for all temperatures, even for the 
more favorable case of 10 X solar metallicity at a C/O ratio of 0.5. A color version of this figure is available in the online journal. 



H 2 photolysis (and other kinetic processes) and due to 
the large overall abundance of CO. The main problem is 
that the reverse of this reaction is also effective at atmo- 
spheric conditions on HD 189733b, due to the large avail- 
able H abundance and the relatively high atmospheric 
temperatures. The forward and reverse reaction rates 
balance each other to a large degree, and although there 
is a small positive net production of CO2 in our dayside 
atmospheric models, the CO2 column abundance never 
achieves the levels needed to exp lain the HST /NICM OS 
secondary-eclipse observations of Sw ain et al.l (|2009al) . It 
is possible that we are missing some key photochemi- 
cal mechanism that effectively converts water to CO2 
on transiting exoplanets, which would b e needed to ex- 
plain a hig h derived CO2/H2O ratio (see lLee et al.ll2012l ; 
iLine et alJ[20TTbh . but given the efficiency of H2O re- 
cycling in a hydrogen-dominated atmosphere (see also 
iMoses et al1l201lHLine et al.ll2010D . it is hard to imagine 
that any chemical mechanism could produce this result 
while hydrogen is still the dominant atmospheric con- 
stituent. 

The problem could potentially b e resolved if the ab- 
sorption attributed to CO2 in the iSwain et al.l (|2009al) 
NICMOS data actually resulted from some other atmo- 
spheric constituent. Species like NH3, HCN, and C2H2 



come readily to mind as candidates, as they are all rel- 
atively abundant in our disequilibrium models, partic- 
ularly for super-solar C/O ratios (for the case of HCN 
and C2H2) or rapid atmospheric mixing at depth (for 
the case of NH3). Hydrogen cyanide has strong bands 
in the ~3-/mi region and much weaker bands in the 
~1.85 /im and ^1.53 fim regions; as such, HCN seems 
unlikely to cause significant absorption in the needed 
1.9-2.2 /im range, despite the relatively large predicted 
column abundance of ~4xl0 19 cm" 2 above 1 bar in 
our nominal C/O = 0.88 disequilibrium model. For the 
same model, NH3 and C2H2 have column abundances of 
~2.5 xlO 19 and ^2.5 xlO 16 cm" 2 , respectively. Acety- 
lene (and HCN) could be more abundant than we derive 
if the upper atmosphere is cooler than we have assumed 
(e.g. JMoses eTaIll201ll ) or if ion chemistry enhances its 
abundance, but the C2H2 bands near ^1.53 and ~2.45 
/jm do not seem good spectral matches to the NICMOS 
data. Ammonia has bands centered bear ~1.97 and 2.2- 
2.4 fxm that could be of interest to the NICMOS data, 
particularly given that models with different assumed 
thermal structures and K zz values (|Moses etall l20lTl) 
predict a significantly larger NH3 abundance than the 
models shown here. Although the NH 3 band at 1.97 
jiva. lines up well with the minimum in flux seen in the 
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NICMOS spectra, the continued absorption in the ~2.1 
/xm region seen in the NICMOS data does not seem con- 
sistent with NH3; however, due to the likely presence of 
NH3 on HD 189733b and other transiting planets, further 
investigation is warranted. At the very least, NH 3 and 
HCN should be considered in both near-IR and thermal- 
IR spectral models for hot Jupiters, as these molecules 
are important disequilibrium products for a variety of 
C/O ratios, as well as important equilibrium products 
for cooler planets and/or for planets with high C/O ra- 
tios. 

Models that contain more ammonia (e.g., because of 
stronger atmospheric mixing or a cooler thermal struc- 
ture at depth) could also help improve the discrepancy 
between our best-fit solutions to the secondary-eclipse 
Spitzer IRAC data as c ompared with the IRS spectr a 
for HD 189733b (see also lMadhusudhan fc Seagenl2009l ). 
A hot lower atmosphere was found necessary to fit the 
high flux in the 3.6 /xm-channel IRAC data i n bot h our 
analysis and that of iMadhusudhan fc Seagerl (|2009[) , but 
that solution also caused excess flux in the 9-12 jum re- 
gion that is i ncons istent with the IRS data (see Fig. 2|). 
iMoses et al.l (|2011l ) demonstrate that the flux in the spec- 
tral models at ~9-12 /xm is considerably reduced when 
opacity from a moderate amount of quen ched NH3 is 
inclu ded in the calculations (see Fig. 13 of IMoses et al.l 
1201 ID. However the very recent Spitzer IRAC analysis 
of lKnutson et all (|2012t ) suggests a very large 7.5a reduc- 
tion in the 3. 6- /mi flux at secondary eclipse, eliminiating 
the need for such a warm lower atmosphere to explain 
the 3.6-/xm data and reducing the model flux in the 9-12 
/zm region to be more in line with the IRS observations 
(see Section [3. 1 .31 and Fig. [6] for more details). 

Although we have focused on secondary-eclipse obser- 
vations in this paper, transit observations can also pro- 
vide compositional information. The C/O ratio is a bulk 
atmospheric property that does not change on observa- 
tional time scales. Atmospheric temperatures, on the 
other hand, change relative quickly through much of the 
photospheric region due to changes in the incidence an- 
gle of the radiation coming from the host star and to 
the relative ly short radiative h eating and cooling time 
scales (e.g... iBarman et aT1l2005t ). Dynamical time scales 
are of order ^10 5 s, which is typically longer than radia- 
tive time scales, allo wing horizontal therm al gradients to 
be maintained (e.g., iShowman et al.ll2009fl . The average 
temperature st ructure therefore cha nges with viewing 
geometry (e.g., iKnutson et al.l [20071. If chemical time 
scales are shorter than dynamical time scales, the com- 
position will also change with viewing geometry. 

Chemical time scales range from very short (<C 10 5 
s) in the deepest, hottest regions below the photosphere 
and in the uppermost regions of the photosphere where 
disequilibrium photochemistry is active, to > 10 5 s in 
the mid-to-lower photosphere. In the mid-to-upper pho- 
tosphere, chemical time scales can either be <10 5 s 
or >10 5 s, depending on stratospheric temperatures. 
When chemical lifetimes are longer than transport time 
scales, which is likely true for a good portion of the 
photosphere on the cooler transiting exoplanets, then 
transport-induced quenching can operate effectively to 
homo genize composition horizon tally as well as vertically 
(e.g.. ICooper fc Showman! 12006) . significantly reducing 



the temperature-related equilibrium compositional vari- 
ations shown in plots like Fig. |T7j However, chemical 
lifetimes will be shorter on warmer planets and for cer- 
tain altitude regions on cooler planets, making compo- 
sitional variations possible between transit observations 
(which are sensitive to the limb atmosphere at the ter- 
minators) and eclipse observations (which are sensitive 
to the fully illuminated dayside atmosphere). Accurate 
predictions of the composition as a function of planetary 
viewing angle will therefore require more sophisticated 
time-variable disequilibrium models, which we delegate 
to future investigations. 

5. THE CASE FOR PLANETARY C/O RATIOS GREATER 
THAN THE HOST STAR 

According to our current understanding of star forma- 
tion, the elemental ratios within a stellar photosphere 
reflect to a large degree the elemental ratios of the pro- 
tostellar nebula in which the star formed, except for a 
slight depletion in heavy elements relative hydrogen due 
to later diffusive set tling in the stellar atmosphere (e.g., 
ILodders et al.ll2009l and references therein). Since plan- 
ets also formed within that same nebula (i.e., protoplan- 
etary disk), one might naively expect giant planets that 
efficiently accrete nebular gas to have elemental ratios 
that match that of the host star and bulk protoplanetary 
disk. In that naive view, giant planets with high C/O ra- 
tios would then be expected to derive from a protoplane- 
tary disk with high bulk C/O ratios, and their host stars 
should be similarly carbon-enriched. Many factors can 
complicate this naive picture (see below), but detailed 
evolution models for WASP-12b suggest that a proto- 
planetary disk with a supersolar C / O ratio is required to 
produc e WASP-12b with C/O > 1 (jMadhusudhan et al.1 
1201111 ). 

Although most stars appear to have near-solar rela- 
tive abundances of carbon and oxygen, there is some 
observational evidence for a C/O-rich tail of the stel- 
lar population distribution, indicating a non-trivial 
percentage of s tars with C /O ratios greater th an 1 
(lEcuvillon et al.l 12004 I2006L iDelgado Mena et all I2010L 
Petigurafc Marcvl 120111: see also iBond et al.l boiolf . 
Gaidosl (|2000l ) points out the C/O ratio of the gas in 
the galactic disk is likely increasing with time, lead- 
ing to a C/O ratio that is larger for younger plane- 
tary systems in our galaxy a nd for those farthe r away 
from the galactic bulge, and iPetigura fc Marcvl ()2011[ ) 
provide observational evidence that stars with plane- 
tary systems have statistically higher C/O ratios than 
those without planets. Both these factors could po- 
tentially explain observatio ns of carbon-rich giant plan- 
ets. On the other hand, iFortnevi ()2012l ) notes that 
stars with C/O ratios greater than 1 have likely been 
overestimated in observational studies (|Ecuvillon et al.| 
200 l|200llDelgado Mena et al.|[20Tol: IPetigura fc Marcvl 
201 If) such that carbon-rich stars (and by correlation, 
carbon-rich protoplanetary disks) are more likely to be 
very rare in the galaxy. In this vein, it is notable that 
the host star WASP-12 appears to have a subsolar C/O 
ratio of 0.40t^7 (|Petigura fc Marcvl 1201 il), in contrast 
to the suggestion of Madhusudh an et all (|2011aD and the 
work presented here that WASP-12b's atmosphere could 
be carbon rich, and in conflict with the naive view of 
similar elemental ratios in a giant planet and its host 
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star. How then could a giant planet develop a C/O ratio 
different from its host star? 

Several processes — some better understood than oth- 
ers — can conspire to produce this result. Disk chemistry 
is not uniform with location and time, and physical and 
chemical fractionations can lead to dif ferent C/O ratios 
in the nebular gas and/or solids (e.g.. iKrot et al. 120001 : 



m the neb ular gas and/or solids (e.g.. IKrot et al. IzUUU 
iCuzzi et a l. 2005; Petaev et al]|2lMlLoddersll2010l) . Ev 



idence from meteorites in our own solar system suggests 
that although the bulk nebula had a solar C/O ratio, 
there were pockets that must have been both oxygen-rich 
and oxygen-poor (at different locations and/or times) 
to account for certain mineral assemblages in different 
meteorite classes, with enstatite chondrites in particu- 
lar suggesting some regions in the solar nebula had C « 
O at the time and loca tion in which the enstatite chon- 
drites were formed (cf. lLarimerl 119751: IKrot et al.1 l200Ct 
iHutson fc Ruzickall2000l : iPasek et al.ll2005h . 

One obvious fractionation process is condensation. 
Condensation of water ice at the "snow line" in a diffu- 
sive or turbulently mixed nebula p rogressively depletes 
water vapor from the inner nebula ([Stevenson fc Lunind 
1988), and this proce ss plus inward drift of ice grains 
from the outer nebula (Stcpinski & Valageas 1997) leads 
to a pile up of water-ice solids near the snow line that 
can greatly facilitate the formation of giant planets (as- 
suming th e core-accretion model for giant-planet fo rma- 
tion; e.g., iMizuno et all 119781 : iPollack et all 119961 ). As 
more water becomes trapped in proto-giant-planet cores 
and other planetesimals, the overall C/O ratio of the 
remaining vapor increases, and water-vapor transport 
to this cold trap c an eventually leave the inner nebula 
deple ted in water (jStevenson fc Lunindll988t iCvr et ail 
19981), although there may be temporary regions of en- 
hanced wa ter vapor where inwardly drifting icy grain s 
evaporate (|Cuzzi fc Zahnldl200l ICiesla fc Cuzzll2006f ). 
The outer nebula will also be depleted in water vapor 
wherever water ice condenses. Several such condensa- 
tion/evaporation fronts for different ices occur at var- 
ious radial distances in the nebula, leading to solids 
with different degrees of carbon and oxygen enrich- 
ments/depletions and vapor with complementary ele- 
mental depletions/enrichments, all of which can influence 
the composition of giant-planet atmospheres (jLodderd 
12001: iDodson-Robinson et al.l l200l lOberg et al.ll20TTI ). 
Solids of sufficient size can decouple from the gas and 
evolve differently in the disk such that although the so- 
lar nebula had a bulk solar C/O ratio initially, the rel- 
ative abundances of carbon and oxygen in the solids 
and gas varies with l ocation and evolves with time (e.g., 
ICiesla fc~C uzzi 2006). The C/O ratio in the atmosphere 
of a giant planet then depends on the location and timing 
of its formation and the subsequent evolutionary history. 

If one assumes that a planet's carbon enrichment is 
largely due to the accretion of solid material, which 
could incidentally also result in planetary metallicities 
enhanced over the host star's, then there is a limited 
region within a protoplanetary disk from which the 
plane t could have accrete d these carbon-enriched solids 
(e.g.. lLodderd [2004. 2010). Water ice condenses beyond 
the H2O snow line (which forms initially near ~5 AU 
in the solar nebula but moves inward with time; see 
ICiesla fc~C uzzi 2006J), whereas carbon can remain pre- 



dominantly as gas - phase CO ou t to its condensation front 
(at > 10-30 AU; iBerginl 12001 IDodson-Robinson et al.l 
2009), resulting in oxygen-enriched solids and carbon- 
enriched gas in the region between the H2O and CO 
condensation fronts. Some additional volatile elements 
can be trap ped in clathrate hydrates within this region 
of the disk (jLunine fc Stevensonlll985l ). but as these ma- 
terials are also water- rich, the solids between the H2O 
and CO condensation fronts are expected to be enriched 
in oxyg en and depleted in ca rbon relative to bulk nebular 
values (jHersant et al.ll2001l ). Therefore, the addition of 
accreted planetesimals from this region of the disk would 
be expected to decrease the C/O ratio of a giant planet's 
atmosphere below the st ellar ratio (e.g .. iGautier et all 
20011: IHersant et al.1 12004). In contrast, ILodderd (|200¥ 
suggests that disks also contain a "tar line" inward of 
the snow line, where refractory organics are stable but 
water ice has evaporated. If a carbon-rich giant planet 
gets its extra carbon from solid material, then it is likely 
from Lodders' tarry solids or from some kind of carbon- 
rich refractory phase that remains in the inner solar sys- 
tem a fter the ice has evaporated (e.g.. lEbel fc Alexanderl 
2011). This scenario suggests that the carbon- rich giant 
planet would have to have accreted these solids during 
or after migration into the inner solar system. 

Note that if the planet were to accrete greatly enhanced 
silicate or metal abundances, a potentially greater frac- 
tion of the oxygen could be tied up in condensed silicates 
and oxides in these planets, leaving the atmosphere above 
these condensate clouds to be oxygen-depleted. However, 
silicate and metal condensates in the protoplanetary disk 
are typically already saturated with oxygen (i.e., metal 
oxides and silicates like MgSiOs and Mg 2 Si04), making- 
it unlikely that the accretion of metals and silicates in 
excess of solar proportions would be a mechanism for 
depleting oxygen in the planets. 

In another scenario, the carbon enrichment in a gi- 
ant planet could result from accreted gas rather than 
from solids. lOberg et alJ (|2011l ) suggest that planets 
like WASP-12b could have formed beyond the H2O snow 
line, with the C/O enhancement simply resulting from 
the enhanced gas-phase nebular C/O ratio that results 
from H2O condensation. For that scenario to work, there 
cannot be much mixing between the protoplanetary core 
(which is presumably water-rich due the planet's forma- 
tion region beyond the snow line) and the atmosphere 
after the planet undergoes its fast runaway gas-accretion 
phase. Moreover, the planet cannot continue to accrete 
large quantities of more oxygen-rich gas as it migrates 
to its present position close to its host star, and it can- 
not accrete a late veneer of oxygen-rich solids that would 
bring its C/O ratio back to the bulk nebular value. If 
the planet becomes massive enough early on to create 
a gap in the disk around its o rbit, then Type II migra- 
tion could ensue (| Ward! [19971) : because the creation of 
the gap can substantial ly slow down the rate at which 
the planet accretes gas (jD'Angelo fc Lubowl 12008). con- 
tinued gas accretion may not be much of a factor for 
these planets once they begin migrating. Accretion of 
solids could still be important, but the current orbital 
positions of the transiting planets lie well within the 
condensation fronts for things like silicate and metals 
within the protoplanetary disk, and evaporation of most 
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condensates might prevent these planets from accreting 
significant quantities of solids once the planets become 
parked at their present orbits. One observational con- 
sequence of this scenario would be a planetary metallic- 
ity the same order as the host star ()Oberg et al.ll20il"l ). 
as superstellar metallicites would imply a larger contri- 
bution from solids (unless photoevaporation of the disk 
is responsible for enhancing virt ually all heavy element s 
in the protoplanetary disk, see iGuillot fc Huesol 12006). 
With this scenario, one would expect C/O ratios near 
unity, as CO is by far the dominant carbon- and oxygen- 
bearing gas in the region between the H2O and CO con- 
densation fronts; note that CO2 seems to be less im- 
ortan t in a column-integrated sense than | Oberget all 
.20111) have assumed (see, e.g.. [Willacy fe Langer 120001 : 
IWillacv fe Wo7^l2l)0a IWalsh et al.ll2012fl . 

There is one other region of the protoplanetary disk 
in which the gas-phase abundance of CO is expected to 
significantly exceed that of water (even assuming bulk 
solar C/O ratios), and that is in the very inner regions 
of the disk (<0.1 AU) where high temperatures (> 1800 
K) begin to affect the stability of water. Water ther- 
mally dissociates at ab out 2500 K, wherea s CO can sur- 
vive to -4000 K (e.g., iNajita et ail 120071 ). Excess oxy- 
gen not tied up in CO would be present as O and + 
in these regions and might be expected to have differ- 
ent vertical distributions than CO due to molecular and 
ambipolar diffusion, with O and + being concentrated 
at higher altitudes in the disk atmosphere than CO. 
At these high altitudes, the O and + might be more 
vulnerable to loss proce sses like stellar- wind stripping 
(|Matsuvama et all I2009D or accelerated accretion onto 
the star ([Elmegreenl 119781) . As such, the gas near the 
mid-plane could be carbon-rich in these warm inner-disk 
regions, with CO the dominant carbon and oxygen con- 
stituent. Disk-chemistry models do not generally extend 
to regions so near the star, but observations of CO over- 
tone emission in young low-mass stellar objects suggests 
that t he CO emission derives from the —0.05 -0.3 AU 
region (|Chandler et al.|[l99H iNaiita et al.ll2000l) . Obser- 
vations of CO and H2O emission obtained for the same 
disk are consistent with the CO being locate d inward of, 
and at higher tempera tures than, the H 2 (|Carr et al.1 
2004; iThi fe Bild 12005). Since hot Jupiters end up in 
these inner-disk regions after whatever migration process 
gets them there, local accretion of gas may also account 
for high atmospheric C/O ratios. Again, since gas accre- 
tion is the culprit for the high C/O ratio in this scenario, 
one might expect the planet's metallicity to be similar 
to the host st ar's, unless H has bee n preferentially lost 
in the nebula (|Guillot fe Hue so 2006). Another potential 
observational test for this scenario would be the N/O ra- 
tio, as N2 would be the only other stable heavy molecule 
abundant under these conditions. 

The fact that we derive C/O ratios near unity for three 
out of four of the hot Jupiters investigated in this paper 
(and the fourth might also be just below 1) seems to 
suggest that the CO-rich gas accretion scenarios from ei- 
ther the hot inner nebula or from beyond the H2O snow 
line (e.g., lOberg et al.l 1201 J ) represent the most likely 
explanation for the atmospheric carbon enrichment, as 
the carbon-rich-solid accretion scenario could potentially 
lead to C/O ^ 1. However, for carbon-rich atmospheres, 



graphite is expected to be a major condensate under 
temperature-pressure conditions in extrasolar giant plan- 
ets, and graphite would thus commandeer the excess 
available car bon, leaving C/O = 1 above the graphite 
clouds (e.g.. lLodders fc Feglevl [T997). The carbon-rich- 
solid scenario is therefore still a viable option, even for 
planets with an apparent C « O. Note, however, that 
WASP-12b is too hot for graphite to form, and if the 
atmosphere had a C/O ratio much greater than 1, we 
would expect to see spectral evidence to that effect. 

Our own giant planets do not shed much light on 
the likelihood of these scenarios. The overall en- 
hanced metallicity for Jupiter and the other giant plan- 
ets suggests that solid planetesimals were instrumen- 
tal in supplying heavy ele ments to our ow n solar- 
system giant planets (e.g., lOwen et"aT1 Il999h . How- 
ever, the ultimate cause of the near-equal enrich- 
ments of most heavy elements measured by the Galileo 
probe is under debate ([Niemann et al.lll998tlWong et al. 
2001 lOwen et al.l 119991: iGautier et al.l 120011: lLodders 



200l lOwen fc Encrenazl 120061: iHersant et all 120041: 



Albert et al.l 120051 : IGuillot fc Huesol 120061 : iMousis et all 
I2009bl ). The fact that the probe entered a dry hot- 
spot region on Jupiter and was not able to mea- 
sure the deep water abundance is unfortunate, as 
the O/H enhancement over solar (or lack thereof) 
could have helped distinguish between different forma- 
tion scenarios (e.g., iLunine et al.l 120041 : lAtreva fc Wongl 
2005). Indirect evidence based on the observed 
CO abundance and disequilibrium-chemistry argu- 
ments is uncertain enough that a range of bulk Jo- 
vian C/0_j^ i y£S_Jspossible from these calcu l ations 
Ccf. iFeglev fc Loddersl Il99l lLodders fc Feglevl l200l 
iBezardet al.l l2002tlVisscher et al.ll2010bl ). although the 
latest models seem to fav or C/O ratios > solar on Jupiter 
(jVisscher fc Mosesll20Tlh . seeming to select against some 
of the high-water-ice and high-clathrate accretion sce- 
narios. Given that Jupiter and transiting planets had 
very different migration histories, it is not clear that their 
carbon and oxygen inventories would have come from a 
similar source. 

In any case, the C/O ratio and bulk metallicity within 
a giant planet's atmosphere provide important clues that 
may help us unravel the puzzle of giant planet formation 
and evolution. Further investigations that can constrain 
these properties from transit and eclipse observations are 
warranted. Ultimately, it will be important to collect a 
statistically significant sample of exoplanet atmospheric 
properties to determine how common carbon-rich giant 
planets are within the exoplanet population. 

6. SUMMARY 

Analyses of transit and eclipse observations sug- 
gest that some transiting extrasolar giant planets have 
unex p ectedly low H 2 Q a bunda n ces (e.g., iSeager et all 
20051: iRichardson et all 120071: iGrillmair et all 1 20071: 
Swain et al.l l2009a|bl: iMadhusudhan fe SeageiT I2OO9I 
2010a| iDesert et all l2009t ISing et all 120091: iLine et al l 
2011bl : lLee et alll2012l : lCowan et alll2012UMadhusudhanl 
20121 ). Atmospheres with solar-like elemental abun- 
dances in thermochemical eq uilibrium are expected t o 
have abundant water (e.g., lLodders fc Feglevl [2002), 
and disequilibrium processes like photochemistry cannot 
seem to deplete water sufficiently in the infrared pho- 
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tospheres of these planets to explain the observation s 
(e.g., see Section 3; lline et aLll2010t iMoses et al.ll201lD . 
While aerosol extinction might help explain the lack of 
water-absorption sign atures in some of the near-infrared 
transit observati o ns (fLecavelier des Etangs et al.l 120081: 
Pont et al.1 120081 : iSing et al.1 120091 1201 It iGibson et all 
2012f) . scattering and absorption from clouds or hazes 



should be less of a factor for the mid-infrared secondary- 
eclipse observations, unless the clouds are optically thick 
in the vertical, contain relatively large particles, and 
are located at pressure l evels within the infrared pho- 
tosphere (e.g.. Hioul 120021 ). Such conditions may not be 
met for the cooler transiting planets like HD 189733b, 
TrES-1, or XO-lb, where major cloud-forming species 
like iron or magnesium silicates will c ondense too deep to 
affect t he se condary-eclipse spectra (jFortnev et al.l[2005l 
l2006bLl2010i ). or for the hottest hot Jupiters like WASP- 
12b, where temperatures remain so hot throughout the 
column that iron and silicates will not condense. Extinc- 
tion from clouds is therefore not expected to be the cul- 
prit for the low derived water abundances from thermal- 
infrared secondary-eclipse observations of many exoplan- 
ets, and if the interpretation is robust, we must look to 
other potential sources for the observed behavior. 

A low water abundance is a natural consequence of 
an atmosphere with a C/O ratio of 1 or greater. Moti- 
vated by the recent derivation of a carbon-rich atmo- 
sphere for WASP- 12b based on seco ndary-eclipse ob- 
servations (jMadhusudhan et all 1201 lal ). we have exam- 
ined the influence of the C/O ratio on the composi- 
tion of extrasolar giant planet atmospheres, both from 
a thermochemical equilibrium standpoint and as a re- 
sult of disequilibrium processes like transport-induced 
quenching and photochemistry. We find that the equilib- 
rium composition of hot Jupiters is very sensitive to the 



C/O ratio (see alsolSeager et al.ll2005t iKuchner fc Seager 
20051 iFortnevI 120051: iLoddersI 201tt IMadhusudhan et al 



2011allbl : IMadhusudhan! I2012t iKopparapu et al.l I2012D 



Carbon monoxide is a major constituent on all hot 
Jupiters whose atmospheres reside on the CO side of 
the CH4-VS.-CO stability field, regardless of the C/O 
ratio (at least for the range 0.1 < C/O < 2), whereas 
H2O and CO2 become significant constituents for C/O 
< 1, and CH4, HCN, and C2H2 become significant con- 
stituents for C/O > 1. Disequilibrium processes do 
not change this conclusion, although photochemistry 
can greatly enhance the HCN and C2H2 abundances, 
and transport-induced quenching can enhance the lower- 
stratospheric abundances of CH4, NH3, a nd HCN, for 
a variety of assumed C/O ra tios (see also IMoses et al.l 
l201lUVisscher fc Mosesli2011f) . Despite rapid photolytic 
destruction by the intense incident UV radiation and de- 
spite kinetic attack by atomic H, water survives in hot- 
Jupiter atmospheres due to the efficient recycling pro- 
cesses that occur in an ^-dominated atmosphere. A low 
derived water abundance for a transiting hot Jupiter may 
therefore be indicative of an atmosphere with a C/O ratio 
greater than solar ( see also Madhusudh an et al.l 12011a: 
lMadhusudhan1l2012f ). 

We compare the results of our thermochemical and 
photochemical kinetics and transport models with 
Spitzer secondary-eclipse data from WASP-12b, XO-lb, 
CoRoT-2b, and HD 189733b — four exoplanets whose at- 
mospheres have been described as having potentially low 



water abundances (e.g., [M adhusud han fc Sea ger 2009; 
IMadhusudhan et all I2011at IMad husudhan! 12012ft . We 
find that disequilibrium models with C/O ~ 1 are con- 
sistent with photometric data from WASP-12b, XO-lb, 
and CoRoT-2b, confirming the possible carbon-rich na- 
ture of these planets. In contrast, spectra from HD 
189733b are more consistent with C/O < 1. In par- 
ticular, our synthetic spectra for HD 189733b compare 
well with Spitzer secondary-eclipse photometric data for 
models with metallicities between 1-5 x solar for mod- 
erately enhanced C/O ratios between 0.88-1.0. These 
fits are not unique: the specific derived C/O ratio that 
provides the best fit to the data will be dependent on 
the assumed thermal structure and atmospheric metal- 
licity, with lower metallicites and/or higher atmospheric 
temperatures allowing lower (i.e., more solar- like) C/O 
ratios to remain consistent with the low inferred water 
absorption from these planets. In those cases, however, 
the relative band ratios in the Spitzer IRAC channels are 
not always well reproduced. 

Our models indicate that hydrogen cyanide is an im- 
portant atmospheric opacity source when C/O > 1, such 
that HCN can even dominate the opacity in the 3.6 
and 8.0 /im Spitzer channels for a variety of conditions. 
For example, we find that HCN is the main absorber in 
the 3.6 and 8.0 /im bands on the very hot WASP- 12b, 
rather than CH4 being resp onsible (as was suggested by 
IMadhusudhan et al.l 1201 lal) or C2H2 being responsible 
(as was suggested by IKopparapu et al.l 12012;): however, 
this result could change for different assumptions about 
the thermal structure and C/O ratio. Acetylene is less 
important overall for the Spitzer bandpasses, but C2H2 
is an important photochemical product in cooler atmo- 
spheres for a variety of C/O ratios and can become a 
significant equilibrium constituent in hotter atmospheres 
when C/O > 1. 

Carbon dioxide is not a dominant constituent in 
our equilibrium or disequilibrium models regardless of 
the C/O ratio (again, for the range 0.1 < C/O < 
2 and near-solar metallicities). A CO2 /H2O ratio 
greate r than unity for HD 18 9733b (see iSwain et all 
| 2009al IMadhusudhan fc Seagerl 12051 ILine et all l20lTbT 
iLee et al would require a metallicity greater than 

2000x solar or a C/O ratio -C 0.1, with both scenarios 
seeming to be inconsistent with the planet's bulk den- 
sity and secondary-eclipse observations. We therefore 
cannot explain the relatively large CO2 a bundance in- 
ferred from fTST/NICMOS observations of (|Swain et all 
l2009al) . In general, we do not expect CO2 to be a major 
constituent in the atmospheres of near-solar-composition 
hot Jupiters. Therefore, spectral modelers may want to 
consider the effects of species like HCN and NH3 instead 
of (or in addition to) CO2 and/or photochemical mod- 
elers will need to come up with new mechanisms to ex- 
plain enhanced CO2 abundances if spectral features on 
hot Jupiters can be tied uniquely to CO2. 

Although an atmospheric C/O ratio > 1 on a giant 
planet could r esult from a carbon-enrich ed protoplane- 
tary disk (e.g.. lMadhusudhan et al.l2011bf ). there are also 
ways in which such a planet could form in a disk that 
possesses more solar-like bulk elemental ratios. Vapor 
condensation within the protoplanetary disk is one sig- 
nificant way in which elements can be fractionated, with 
the solids and vapor ending up with very different C/O 
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ratios. Carbon-rich giant planets could have accreted 
water-ice-poor, carbon-rich solids during their migration 
through the disk regions inward of th e snow line (e.g., 
lLoddersll2004tlEbel fc Alexanderll2011| ). or they could get 
their carbon enrichment from the accretion of CO-rich, 
H^O-poor gas, either during the runaway accretion phase 
as they formed in the region between the water and CO 
cond ensation fronts in the outer disk (e.g., IQberg et al.1 
2011) or in situ from the inner disk once their migra- 
tion stopped very close to the host star. Regardless of 
whether these specific scenarios are viable or not, the 
atmospheric C/O ratio and bulk metallicity provide im- 



portant clues regarding the formation and evolution of 
giant planets, and it is hoped that future spectral obser- 
vations can reveal how common carbon-rich giant planets 
are within the exoplanet population. 
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